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Introduction 

^Htork  on  this  program  during  the  past  year  has  been  directed  towards: 

(1)  analysis  of  surface  damage  that  occurs  during  sliding  contact;  (2) 
measurement  at  elevated  temperatures  of  friction  between  ceramics;  (3) 
investigation  of  point  and  line  damage  that  occurs  during  sharp  contact. 

Surface  Damage  Studies 

The  surface  damage  studies  were  conducted  on  the  same  equipment  that  was 
developed  on  this  program  to  study  strength  degradation.  The  equipment  is 
designed  to  measure  the  sliding  coefficient  of  friction  of  a  sphere  against  a 
flat  specimen,  and  is  capable  of  operating  at  temperatures  as  high  as  800°C. 
With  modification,  the  equipment  Is  capable  of  operating  at  temperature  up  to 
1500°C. 

Initial  studies  were  conducted  on  soda-lime  silica  glass.  Examination  of 
the  wear  tracks  Indicates  that  plastic  flow  occurs  at  a'l  temperatures  when 
glass  slides  upon  glass.  As  the  temperature  Is  Increased  the  zone  of  plastic 
flow  and  the  friction  are  observed  to  increase,  primarily  as  a  result  of 
plowing  of  the  glass  sphere  through  the  glass  specimen.  At  elevated 
temperatures  crack  formation  Is  most  severe  at  the  point  where  sliding 
initiates.  The  size  of  the  crack  can  be  correlated  with  the  higher  frictional 
forces  required  to  Initiate  sliding  at  elevated  temperatures.  At  the  fictive 
temperature  large  cracks  are  observed  at  the  start  of  the  wear  track;  however, 
little  or  no  cracking  Is  observed  along  the  wear  track.  By  contrast,  crack 
formation  is  much  more  uniform  along  the  wear  track  at  room  temperature.  The 
relevance  of  the  crack  formation  to  the  strength  of  the  glass  is  currently 
being  studied. 
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Friction  Data 


Friction  data  were  obtained  from  a  plate  on  plate  apparatus  that  was 
built  to  measure  friction  where  mechanical  damage  does  not  occur  during 
sliding.  Results  from  the  apparatus  indicate  that  the  coefficient  of  friction 
between  two  solids  increases  dramatically  as  the  temperature  is  increased  to 
levels  at  which  molecular  transport  can  occur  easily.  For  soda-lime  silica 
glass  and  a  70  percent  lead  silica  glass,  the  coefficient  of  friction  was  ^.3 
from  room  temperature  to  -v200oC  below  the  fictive  temperature.  As  the 
temperature  increased  the  coefficient  of  friction  increased  to  -v 1  at  the 
fictive  temperature  of  the  glass  and  then  decreased  again  at  higher 
temperatures.  Studies  on  silica  glass,  aluminum  oxide  and  silicon  nitride  (NC 
132)  gave  no  indication  of  such  dramatic  behavior  because  of  the  higher 
temperatures  needed  for  molecular  transport.  This  part  of  the  work  was 
conducted  by  John  S.  Nadeau,  who  was  partially  supported  by  this  contract 
during  the  past  year. 

Studies  using  glass  spheres  on  soda-lime  silica  glass  gave  indications  of 
a  similar  increase  in  the  coefficient  of  friction.  This  increase  could  not  be 
attributed  to  adhesion  effects  since  measurements  using  normally  loaded  glass 
spheres  against  flat  plates  gave  no  indication  of  adhesion  until  the 
temperature  exceeded  the  fictive  temperature  of  the  glass.  Apparently, 
tangential  as  well  as  normal  forces  are  needed  for  the  observed  frictional 
effects. 

Point  and  Line  Damage 

Studies  of  rate  effects  on  hardness  of  materials  were  conducted  using  a 
dynamic  hardness  machine  that  was  constructed  in  an  earlier  part  of  the  ONR 
program.  Tests  were  applied  to  single  crystal  copper,  to  polycrystalline 


2 


tungsten  and  to  soda  lime  silica  glass.  The  contact  times  ranged  from  ^10  ms 

3 

to  10  s.  A  distinct  decrease  in  hardness  was  observed  over  this  range: 

12.5%  for  copper;  13.3%  for  glass;  21.3%  for  tungsten.  The  results  of  this 
study  have  been  accepted  for  publication  in  the  Journal  of  Materials  Science 
(Letters). 

In  addition  to  the  above,  a  theory  was  formulated  for  the  general  fatigue 
response  of  flaws  in  brittle  materials  to  residual  contact  stresses  associated 
with  both  point  and  line  indentations.  These  residual  stresses  can  play  a 
dominant  role  in  the  ensuing  failure  mechanics.  Analytical  solutions  were 
obtained  for  the  specific  case  of  static  fatigue.  The  resulting  relation 
between  lifetime  and  failure  stress  has  the  same  power-law  form  as  the 
conventional  solution  for  residual  stress  free  flaws.  However,  while  the  form 
of  the  lifetime  relation  remains  invariant,  the  values  of  the  exponents  and 
coefficients  depend  on  the  type  of  flaw  (line  versus  point  flaw)  and  on  the 
presence  of  residual  stresses.  Explicit  conversion  formulae  are  given  to 
transform  "apparent"  velocity  parameters  into  "true"  parameters.  This  work 
has  been  submitted  for  publication  in  the  Journal  of  the  American  Ceramic 
Society. 
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Effect  of  material  parameters  on  the  erosion 
resistance  of  brittle  materials 
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Erosion  data  are  compared  with  two  theories  that  have  been  suggested  to  explain  the 
erosive  behaviour  of  solids.  A  dimensional  analysis  is  applied  to  the  variables  that  are 
important  to  erosion,  and  a  multivariate,  linear  regression  analysis  is  used  to  fit  the  data 
to  the  dimensional  analysis.  The  results  of  the  linear  regression  analyses  are  compared 
with  the  two  theories  in  order  to  evaluate  the  applicability  of  these  theories  to  erosion. 
Although  semi-quantitative  agreement  of  the  data  with  the  theories  is  obtained,  some 
discrepancies  are  apparent.  In  particular,  the  dependence  of  erosion  rate  on  hardness  and 
critical  stress  intensity  factor  is  greater  than  predicted  by  either  of  the  two  theories. 
These  discrepancies  are  attributed  primarily  to  microstructural  aspects  of  erosion  that  are 
not  modelled  by  either  of  the  theories. 


1.  Introduction 

Erosion  of  brittle  materials  by  hard,  solid  particles 
is  a  complex  process  in  which  material  is  lost  from 
the  target  surface  by  brittle  fracture  [1,  2],  The 
sizes  and  types  of  cracks  that  form  in  the  target 
surface  during  impact  have  been  studied  extensively 
and  have  been  shown  to  depend  on  several  factors: 
these  include  particle  shape,  mass  and  velocity, 
and  target  material  hardness  and  toughness.  At  low 
velocities  well-developed  crack  systems  form 
at  the  impact  site:  cone  cracks  are  formed  by 
rounded  (“blunt”)  particles  [3-5],  lateral  and 
medium  cracks  by  angular  (“sharp”)  particles  [3 , 
6-8].  At  extremely  high  velocities  the  appearance 
of  the  impact  surface  is  affected  by  the  ejection  of 
material  from  the  target  surface  as  the  particle 
plows  into  the  surface,  and  by  severe  cracking  and 
chipping  of  the  surface  after  the  particle  has  left 
the  impact  site  [1,  4,  9,  10].  The  types  of  crack 
systems  that  are  formed  during  impact  and  the 
conditions  that  control  their  formation  have  been 
discussed  extensively  in  the  references  cited 
above. 

Plastic  deformation  also  plays  an  important 
role  in  the  erosion  process.  Thus,  a  detailed 
examination  of  impact  sites  in  brittle  materials 
indicates  that  a  zone  of  intense  plastic  defor¬ 


mation  forms  during  contact,  beneath  the 
immediate  area  of  the  contact  (7,  11-13], 
Residual  stresses  associated  with  the  plastic 
zone  force  small  cracks,  known  as  lateral  cracks,  to 
grow  from  the  impact  site.  Initially,  these  cracks 
grow  parallel  to  the  target  surface,  but  then  curve 
towards  and  eventually  intersect  with  the  surface 
resulting  in  a  loss  of  material  from  the  target. 
Because  of  this  behaviour,  the  erosion  process  in 
brittle  materials  is  viewed  by  many  investigators  as 
an  elastic-plastic  event,  the  plastic  deformation  at 
the  impact  site  being  the  prime  driving  force  for 
the  surface  fracture  that  results  in  material  loss 
during  erosion  ]1,  2,  10].  While  this  view  of 
erosion  may  be  over  simplified  considering  the 
complexity  of  the  process,  it  has  been  used  to 
model  the  erosion  process  and  to  develop  equations 
that  predict  erosion  rates  as  a  function  of  projectile 
and  target  parameters  that  are  known  to  influence 
the  erosion  process  is  brittle  materials. 

Two  elastic -plastic  theories  have  been  developed 
to  explain  the  erosion  of  brittle  solids.  Both  are 
based  on  the  assumption  that  lateral  cracks  grow 
in  a  quasi-static  manner  as  a  result  of  residual 
stresses  introduced  by  the  impact  event.  In  both 
theories,  the  size  of  the  lateral  cracks,  c,  are 
assumed  to  be  determined  by  the  following 
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relation  [14]*: 

P/cin  ~fiKc  (1) 

where  P  is  the  maximum  normal  load  during 
impact,  Kc  is  the  critical  stress  intensity  factor, 
and  0  is  a  nondime nsional  constant.  The  volume 
of  material  removed  during  erosion,  V,  is  deter¬ 
mined  from  the  size  of  the  lateral  crack,  c,  and  the 
depth  of  the  crack,  d,  beneath  the  target  surface 

V  =  itc3d.  (2) 

Since  the  impact  sites  are  assumed  to  be  non¬ 
interacting,  the  total  wear  volume,  W,  is  just 
the  summation  of  the  volumes  resulting  from 
the  individual  impact  events. 

The  two  elastic -plastic  wear  theories  differ  in 
their  assumed  dependence  of  impact  load,  P, 
on  the  kinetic  and  material  parameters  that  are 
important  to  erosion.  The  theory  developed  by 
Evans  et  al.  [10]  includes  dynamic  stress  wave 
effects  in  the  calculation  of  P.  A  spherical  particle 
is  assumed  to  penetrate  into  a  target  without 
distortion;  the  contact  pressure  is  assumed  to  be 
equal  to  the  dynamic  pressure  that  occurs  when 
the  particle  first  hits  the  target  surface.  The  depth 
of  penetration  is  determined  from  the  time  of 
contact  and  the  mean  interface  velocity,  both 
of  which  are  calculated  from  a  one-dimensional 
analogue.  The  final  expression  for  the  erosion 
rate,  W,  is 

x  [(ZtZp),/3/(Zt1,J  +zp''T3]  O) 
where  v0  is  the  initial  particle  velocity,  R  and  p 
are  the  particle  radius  and  density,  respectively, 
Ke  and  H  are  the  target  toughness  and  hardness, 
respectively,  and  Zt  and  Zp  are  the  impedances  for 
the  target  and  the  particle,  respectively,  H  is  the 
hardness  and  V0  is  the  initial  particle  velocity.  The 
term  within  the  brackets  varies  by  less  than  10% 
for  the  materials  used  in  the. current  study,  and 
therefore  will  be  considered  to  be  a  constant  for 
the  purpose  of  this  paper.  Hence,  the  equation  for 
the  erosion  rate  reduces  to 

w<x  wS.*flS.7pl.S/r1.2S  (4) 

A  quasi-static  formulation  of  the  erosion 


problem  is  based  on  work  by  Wiederhorn  and 
Lawn  [15],  in  which  the  kinetic  energy  of  the 
particle  is  assumed  to  be  absorbed  completely 
by  plastic  flow  when  a  particle  impacts  the  sur¬ 
face.  From  this  assumption,  both  the  maximum 
force  during  contact  and  the  maximum  depth  of 
penetration  can  be  calculated.  Assuming  that  the 
lateral  cracks  generate  at  a  distance  beneath  the 
surface  that  is  equal  to  the  maximum  depth  of 
particle  penetration,  the  following  equation  for 
the  erosion  rate  is  derived  [2] 

W<*vl*R3V*K;l-3H°"  (5) 

The  forms  of  the  two  erosion  theories  presented 
above  are  similar  in  that  they  express  the  erosion 
process  by  a  power  law  dependence  of  erosion 
rate  on  both  particle  (v0,  R,  p)  and  target  (Ke,H) 
properties.  Although  the  same  properties  are  used 
in  both  theories,  the  exponents  for  velocity, 
particle  density,  and  hardness  differ.  A  comparison 
of  these  theories  with  experimental  results  on 
erosion  indicates  that  the  theories  are  reasonably 
consistent  with  experiment  with  regard  to  the 
exponents  for  velocity  and  particle  size  [2], 
The  effect  of  particle  density  on  erosion  has 
not  been  investigated  in  any  systematic  manner, 
so  that  there  is  no  way  of  knowing  if  the  exponents 
given  in  Equations  4  and  5  are  correct.  A  study  of 
the  effect  of  hardness,  H,  and  fracture  toughness, 
Kc,  on  erosion  has  recently  been  conducted  on  a 
series  of  ceramics  by  Evans  et  al.  [10]  and  by 
Gulden  [16].  The  data  obtained  by  Evans  et  al. 
[10]  suggest  a  greater  dependence  of  erosion 
rate  on  Ke  and  H  than  is  predicted  by  their  theory. 
Aside  from  these  studies,  however,  there  have  been 
no  systematic  investigations  of  the  effect  of  Kc 
and  H  on  the  erosion  rate  of  brittle  materials. 

In  this  paper,  the  erosion  of  dense  brittle 
materials  is  studied  in  order  to  assess  the  validity 
of  the  erosion  theories  represented  by  Equations 
4  and  5.  Of  particular  interest  to  this  study  are 
the  particle  velocity  and  the  material  parameters 
Ke  and  H.  The  results  of  our  study  will  show  that 
while  both  theories  provide  a  qualitative  descrip¬ 
tion  of  the  erosion  data,  neither  theory  is  quan¬ 
titatively  correct .  The  reason  for  these  differences 
seems  to  lie  in  the  simplifying  assumptions  made 


‘This  relation  is  concerned  with  the  formation  of  radial  cracks  from  a  sharp  indentation.  For  the  relation  to  be  valid, 
the  crack  must  be  large  relative  to  the  size  of  the  indentation.  The  use  of  this  relation  to  describe  lateral  crack  formation 
is  based  on  work  by  Evans  et  al.  1 101  who  showed  experimentally  that  the  size  of  the  two  cracks  were  proportional. 
Hence,  the  use  of  Equation  1  to  describe  erosion  phenomenon  has  its  basis  in  empirical  investigations,  but  has  no 
theoretical  justification. 
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TABLE  1  Properties  of  target  materials  used  in  erosion  study 


Material 

Young’s 

modulus 

(GPa) 

Hardness 

(GPa) 

Toughness,  Kc 
(MPam1'1) 

Mierostructure 

Hot-pressed 
silicon  nitride 

317  ( 17 1 

19.9  [18] 

5.0  [19] 

Fully  dense  ~  1  jum 
grain  size 

Hot-pressed 
silicon  carbide 

466  {201 

29.4  [18] 

4.0(21] 

Fully  dense  ~  1  to  2  Mm 
grain  size 

Hot-pressed 
aluminium  oxide 

425  [22] 

22.0  [18] 

4.0  [23] 

Fully  dense  -  3  to  4  Mm 
grain  size 

Sintered 

aluminium  oxide 

425  [22) 

21.7  [18] 

2.2  [24] 

Fully  dense  30 Mm 
grain  size 

Sapphire 

425  [22] 

21.7  [18] 

2.2  [24] 

Single  crystal  {lOTl} 
plane 

Silicon 

168 [25] 

10.6  [181 

0.7  (26) 

Single  crystal  {l  1 2} 
plane 

Silica  glass 

75  [27] 

8.7  [18] 

0.77  [28] 

C7940 

Soda-lime - 
silica  glass 

75  [27] 

6.3  [18] 

0.75  [28] 

C0800 

Magnesium 

oxide 

330  [22] 

8.0  [29] 

2.6  [18] 

Fully  dense  ~  10  to  15  Mm 
grain  size 

in  both  theories  of  erosion.  As  will  become 
apparent,  details  of  the  microstructure  and  material 
interaction  during  impact  affect  erosion  in  ways 
not  fully  accounted  for  by  the  present  models 
of  erosion. 

2.  Experimental  procedure 

The  target  materials  used  in  the  present  investi¬ 
gation  provided  a  reasonably  wide  range  of  target 
properties  (Kc  and  H)  and  microstructure  for 
study  (Table  I).  Examination  of  the  impact  area 
by  transmission  electron  microscopy  showed  that 
all  of  the  materials  selected  for  investigation 
deformed  plastically  when  subjected  to  impact 
[13]*.  The  cracks  that  were  generated  by  the 
impacting  particles,  while  originating  from  within 
the  deformed  zone,  exhibited  no  evidence  of 
localized  plastic  deformation  at  the  crack  tip, 
and  accordingly  propagated  in  a  brittle  mannert. 
Therefore,  these  materials  fit  within  the  framework 
of  the  theories  discussed  above . 

The  particles  used  for  erosion  measurements 
were  1  SO^tm  SiC  abrasive  grains.  Because  of  the 
hardness  of  these  particles  it  was  felt  that  they 


simulated  the  hard,  non-yielding  particles  assumed 
to  be  responsible  for  erosion  in  the  theories  used 
to  derive  Equations  4  and  5.  To  achieve  a  uniform 
particle  size  for  investigation,  all  particles  were 
sieved  between  an  80  and  120  mesh  screen  before 
being  used  in  erosion  investigations.  The  particles 
that  were  used  passed  through  the  80  mesh  screen, 
but  were  retained  by  the  120  mesh  screen. 

The  erosion  apparatus  used  in  this  study  has 
been  described  previously  [30].  Briefly,  the 
equipment  was  desgined  to  feed  abrasive  particles 
into  a  high  velocity  air  stream,  which  propelled 
the  particles  against  the  specimen  surface  (Fig.  1). 
The  particles  were  accelerated  by  passing  them 
through  a  tungsten  carbide  nozzle  ~5cm  long 
and  0.16  cm  in  internal  diameter.  The  acceler¬ 
ation  of  the  particles  to  high  velocity  is 
accomplished  within  the  nozzle.  The  particle— air 
mixture  is  passed  through  a  ceramic  tube, 
2  cm  in  diameter,  to  obtain  a  relatively  uniform 
beam  of  abrasive  particles.  High  temperatures 
can  be  achieved  by  feeding  a  propane— oxygen 
mixture  through  a  ring-burner  into  the  top  of  the 
ceramic  tube.  The  high  velocity  particle-air 


*The  technique  of  transmission  electron  microscopy  was  applied  only  to  the  crystalline  materials  used  in  this  investi¬ 
gation.  With  regard  to  their  deformation  and  fracture  properties,  however,  other  studies  suggest  that  the  behaviour  of 
glass  is  similar  to  that  of  crystalline  materials  (see  Lawn  et  at.  [  1 3  ]  for  a  discussion  of  this  point). 

'With  the  exception  of  MgO  for  which  dislocations  can  probably  be  generated  at  crack  tips. 
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Figure  1  Schematic  diagram  of  erosion  equipment  (after 
Wioderhom  and  Roberts  [30) ). 


stream  sucks  the  flame  from  the  burner  into  the 
ceramic  tube  to  produce  temperatures  as  high  as 
1200°  C. 

The  particle  velocities  were  measured  by  using 
a  time-of-flight  technique  developed  originally 
by  Ruff  and  Ives  [31].  In  this  technique,  two 
discs  rotate  on  a  common  axis  which  is  parallel 
to  the  direction  of  the  erosive  gas  stream.  The 
disc  closest  to  the  exit  port  of  the  erosion 
apparatus  contains  a  slit,  which  permits  particles 
to  past  through  the  disc  and  impinge  on  the 
second  disc.  For  a  fixed  rotational  speed,  the 
position  of  the  erosion  mark  on  the  second  disc, 
relative  to  the  position  of  the  slit  on  the  first 
disc,  establishes  the  particle  velocity. 

The  sensitivity  of  the  double  disc  technique 
to  measure  partide  velocities  was  improved 
during  the  course  of  the  present  study  by  mount¬ 
ing  partially  silvered  glass  microscope  slides  on  the 
second  disc  directly  below  the  slit  on  the  first  disc 
[32].  Particles  that  impinged  on  the  glass  slide 
formed  impact  damage  that  destroyed  the  reflec¬ 
tivity  of  the  microscope  slide  in  the  immediate 
area  of  impact.  As  a  consequence,  the  position 
of  the  erosion  marks  on  the  microscope  slides 
were  easily  observed  and  measured.  Quantitative 
optical  microscopy  was  used  to  improve  the 
accuracy  of  measuring  the  position  of  the  erosion 
marks  on  the  microscope  slides.  The  original 
technique  was  further  improved  by  first  rotating 
the  disc  in  one  direction  to  obtain  an  erosion 


mark  and  then  in  the  opposite  direction  to  obtain 
a  second  erosion  mark.  This  procedure  doubled 
the  distance  between  marks,  thus  improving  the 
accuracy  of  the  velocity  measurements. 

The  specimens  used  in  this  study  were  ~  1 .25 
cm  square  plates  approximately  0.6cm  thick. 
They  were  mounted  on  a  support  arm  and  held 
with  their  wide  face  normal  to  the  stream  of 
erosive  particles.  Specimens  were  exposed  to 
a  fixed  mass  of  erosion  particles,  which  ranged 
from  25  to  400  g  depending  on  the  target  material 
and  the  particle  velocity  selected  for  study.  The 
mass  lost  by  the  target  during  each  experiment 
was  measured  to  at  least  1%  accuracy  using  an 
analytical  balance.  The  erosion  rate  was  calcu¬ 
lated  from  the  fraction  of  particles  that  inter¬ 
sected  the  specimen.  The  number  of  particles 
impacting  the  target  was  estimated  from  the 
mass  of  abrasive  used  and  the  mean  particle 
size  of  the  abrasive  (approximating  the  particles 
as  spheres).  Finally,  the  volume  loss  per  particle 
impact  (i.e.  the  erosion  rate)  was  calculated  from 
the  mass  lost  from  the  specimen  per  particle 
impact ,  and  the  target  density . 

3.  Results 

The  results  of  our  studies  are  shown  in  Figs.  2  to 
4*.  In  each  case  the  log  of  the  erosion  rate 
(expressed  as  volume  lost  per  particle  impact)  is 
plotted  against  the  log  of  the  particle  velocity  .  Fig. 
2  presents  the  results  obtained  at  room  tempera¬ 
ture.  The  erosion  rate  was  measured  for  velocities 
ranging  from  37  to  94  m  sec'1 ,  for  the  nine  target 
materials  used  in  the  present  study.  The  erosion 
data  shown  in  Fig.  2  fit  a  power  law  function  as 
expressed  by  Equations  4  and  5.  The  slopes  of  the 
curves  at  room  temperature  ranged  from  1.9  for 
hot-pressed  silicon  carbide  to  2.9  for  silicon  and 
silica  glass.  The  standard  error  of  the  slopes  ranged 
from  ~  0.003  to  ~0.25  with  a  mean  value  of 
~  0.1  (Table  II),  which  indicates  that  at  the  95% 
confidence  level  and  for  two  degrees  of  freedom, 
a  difference  in  slope  of  ~0.4  is  significant.  With 
the  exception  of  the  hot-pressed  silicon  nitride, 
the  values  of  these  slopes  are  similar  to  those 
reported  by  other  investigators  on  similar  materials 
[33—38].  The  slope  of  the  hot -pressed  silicon 
nitride  was  about  one -half  that  reported  earlier 
by  Gulden  [16],  The  erosion  rate  of  the  target 
materials  shown  in  Fig.  2  decreases  as  the  tough- 


The  data  used  in  these  figures  are  summarized  in  Ap  'ndix  B. 


Figure  2  Erosion  of  brittle  materials  at  25°  C,  normal 
incidence  impact,  150  pm  SiC  particles.  For  clarity  the 
errors  given  in  Table  IB  have  been  left  off  the  figure. 


ness  of  the  target  material  increases,  a  finding 
that  provides  qualitative  support  for  the  erosion 
theories  described  by  Equations  4  and  5.  A  quan¬ 
titative  comparison  of  the  two  theories  with  the 
data  will  be  made  in  a  later  section  of  this  paper. 

The  '  ive  wear  data  collected  at  500  and 
1000°  C  are  shown  in  Figs.  3  and  4,  respectively, 
for  several  of  the  target  materials  used  in  the 
present  study.  The  data  shown  in  these  figures 
are  similar  to  those  obtained  at  room  temper¬ 
ature.  However,  for  some  of  the  materials  the 
slopes  of  the  curves  at  elevated  temperatures 
were  significantly  greater  than  those  obtained 
at  room  temperature.  Although  the  relative 
position  of  the  erosion  curves  on  the  graph  was 
roughly  the  same  at  elevated  temperature  and 
room  temperature,  small  systematic  differences 
in  erosion  behaviour  were  obtained  for  some  of 
the  materials.  Thus, elevated  temperatures  appeared 
to  slightly  enhance  the  erosion  rate  of  silicon, 
and  hot-pressed  silicon  nitride  at  the  higher 
velocities,  whereas  the  erosion  rate  of  glass, 
sapphire,  and  sintered  aluminium  oxide  was 
reduced  at  the  lower  velocities.  The  results  of 
the  present  study  were  similar  to  those  reported 
earlier  by  the  present  authors  on  a  smaller  set  of 
data  [7], 


Particle  velocity  (msec'1) 

Figure  3  Erosion  of  brittle  materials  at  500°  C,  normal 
incidence  impact,  150  urn  SiC  particles.  For  clarity  the 
errors  given  in  Table  IB  have  been  left  off  the  figure. 

4.  Discussion 

The  purpose  of  this  paper  is  to  present  data  that 
can  be  compared  with  the  elastic-plastic  theories 
(Equations  4  and  5)  that  have  been  developed  to 
explain  the  erosion  of  brittle  materials.  In  particu¬ 
lar  the  erosion  data  were  used  to  evaluate  the 
exponents  of  v0,  Kc,  and  H,  which  were  then  com¬ 
pared  with  those  given  in  Equations  4  and  5 .  Since 
Kc  and  H  are  determined  by  the  properties  of  the 
target  material,  they  cannot  be  varied  indepen¬ 
dently,  and  hence  have  to  be  compared  with  the 
theories  in  combined  form:  K~l3H~ 0-25  for 
Equation  4;  for  Equation  5.  A  second 

way  of  comparing  the  exponents  of  Equations 
4  and  5  with  the  experimental  data  is  by  first 
expressing  these  two  equations  in  dimensionless 
form  through  the  use  of  a  dimensional  analysis 
and  then  fitting  the  dimensionless  equation  to 
the  experimental  data  to  obtain  the  exponents. 
Both  of  these  techniques  will  be  used  in  this 
paper. 

4.1 .  Velocity  exponents 
The  velocity  exponents  obtained  in  this  paper 
are  summarized  in  Table  II,  and  can  be  com¬ 
pared  with  other  data  reported  in  the  literature 
for  similar  materials  (Table  III).  As  can  be  seen 
from  these  tables,  data  obtained  in  the  present 


Figure  4  Erosion  of  brittle  material!  at  1000s  C,  normal 
incidence  impact,  150  pm  SiC  particles.  For  clarity  the 
errors  given  in  Table  IB  have  been  left  off  the  figure. 

study  are  reasonably  consistent  with  those  reported 
by  other  authors.  Most  of  the  differences  between 
the  results  shown  in  Tables  II  and  III  are  believed 
to  be  due  to  small,  systematic,  interlaboratory 
differences  in  experimental  technique.  The  results 
on  hot-pressed  silicon  nitride,  however,  differ 
significantly  from  our  own  because  of  the  large 


difference  in  velocity  exponent  (4  versus  2.2) 
obtained  in  the  two  studies. 

The  velocity  exponents  from  Table  II  cluster 
more  closely  about  the  value  (2.4)  predicted  by 
the  quasi-static  model  of  erosion,  Equation  5,  than 
the  value  (3.2)  predicted  by  the  dynamic  model  of 
erosion,  Equation  4.  This  conclusion  has  to  be 
tempered  by  the  fact  that  the  velocity  exponent  of 
V  increases  as  the  temperature  is  increased  and 
tends  to  fall  between  the  two  predicted  values. 
Furthermore,  recent  studies  on  silicon  by  Scatter- 
good  and  Routbort  [36]  suggest  that  the  velocity 
exponent  increases  as  the  particle  size  decreases. 
Hence,  for  small  particles,  the  trend  is  toward 
better  agreement  between  the  dynamic  theory  of 
erosion  and  experimental  measurement.  These 
dependences  of  velocity  exponent  on  tempera¬ 
ture  and  on  particle  size  are  not  predicted  by 
either  theory. 

4.2.  Target  parameters 
The  erosion  data  presented  in  Fig.  2  are  compared 
with  the  material  parameters  H  and  Ke  in  Fig.  S. 
Fig.  Sa  compares  the  erosion  data  with  the  dynamic 
theory  of  erosion;  Fig.  Sb  compares  the  erosion 
data  with  the  quasi-static  theory  of  erosion.  With 
the  exception  of  MgO,  the  data  on  both  figures 
plot  as  straight  lines,  lending  credence  to  the 
suggested  theories  of  erosion.  However,  both  sets 
of  data  are  represented  by  lines  with  slope:  greater 
than  1,  the  theoretically  expected  slope.  The 
empirical  slope  for  the  dynamic  erosion  theory, 
~  1 .2,  is  closer  to  the  expected  slope  of  1  than  is 
the  slope,  ~1.5,  for  the  quasi-static  theory  of 


TABLE  II  Velocity  exponent!  for  eroiion  data:  normal  incidence 


Material 

25°  C 

Temperature 

500°  C 

1000°  C 

Magnesium  oxide,  polycrystalline 

2.2 

- 

- 

Soda-Ume -silica  glass 

2.5  (0.12)* 

3.5  (0.20) 

- 

Vitreous  silica 

2.9 

3.0 

- 

Sapphire 

2.3(0.10) 

2.4  (0.25) 

3.3  (0.03) 

Sintered  aluminium  oxide,  30  pm 

2.3  (0.003) 

2.8  (0.09) 

2.7  (0.15) 

Hot-pressed  aluminium  oxide,  3  to  4  run 

2.3  (0.03) 

2.1  (0.04) 

2.3  (0.11) 

Silicon 

2.9  (0.03) 

3.8 

3.4 

Hot-pressed  silicon  carbide 

1.8(0.16) 

- 

- 

Hot-pressed  siUcon  nitride 

2.1  (0.08) 

2.5  (0.03) 

2.4  (0.20) 

*The  number!  in  parentheaei  give  the  standard  error  for  the  value  of  the  velocity  exponent,  which  was  determined 
by  a  Unear  regression  analysis  of  the  mean  wear  values  given  in  Table  IB.  For  exponents  that  were  determined  from 
only  two  wear  values,  no  standard  error  is  given. 
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TAB  L  E  III  Velocity  exponents  for  erosion  data:  normal  incidence 


Target  material 

Erosion  particles 

Exponent 

Reference 

Soda-lime -silica  glass 

SiC,  120  grit 

3.0 

133) 

MgO  (96.5%) 

SiC,  120  grit 

2.7 

[331 

Al,Oj  (99.5%) 

SiC,  120  grit 

2.6 

[33] 

Pyrex  glass 

AljO,  30 

2.2 

[34] 

10  gm 

2.7 

Hot-pressed  SijN, 

SiC  8  to  940  gm 

4.0 

[35 ) 

Silicon 

A12Oj  23 

3.4  to  2.6 

[36] 

to  270  um 

depending 
on  particle 
size 

Reaction  bonded  SiC 

AljO,  130Mm 

2.3 

[37] 

270  Mm 

2.0 

Hot-pressed  SiC 

AljO,  130 Mm 

1.8 

1 38] 

270  pm 

1.5 

erosion.  Similar  values  of  slopes  were  observed 
by  Gulden  [16]  in  an  erosion  study  on  a  different 
set  of  materials.  Hence,  from  this  type  of  analysis, 
the  dynamic  theory  of  erosion  appears  to  provide 
a  somewhat  better  fit  to  the  erosion  data  than 
does  the  quasi-static  theory. 

The  data  for  MgO  are  not  consistent  with  the 
data  obtained  for  the  other  materials,  undoubtedly 
because  of  the  type  of  impact  damage  formed  in 
the  surface  of  this  material.  The  MgO  cracked 
along  the  grain  boundaries  in  the  vicinity  of  the 
impact  site,  so  that  each  impact  event  formed  a 
loosely  connected  aggregate  of  grains  that 


surrounded  the  impact  site.  These  damaged  areas 
were  easily  removed  from  the  surface  during 
erosion,  resulting  in  a  higher  rate  of  erosion  than 
predicted  theoretically.  In  essence, the  mechanism 
of  erosion  for  MgO  differed  markedly  from  that 
for  other  materials.  Erosion  of  polycrystalline 
MgO  probably  does  not  fit  the  lateral  chipping 
models. 

4.3.  Dimensional  analysis 
Dimensional  analysis  [39]  is  an  alternative  method 
of  obtaining  relationships  between  the  parameters 
that  affect  erosion.  While  not  providing  a  specific 


(a)  "°25  lb)  /rc-1.3J*0.11 

Figure  5  Comparison  of  the  erosion  data  with  the  theories  given  by  (a)  Equation  4  (10)  and  (b)  Equation  5(2).  The 
volume  lost  per  particle  impact  was  selected  from  Fig.  2  at  a  velocity  of  63  m  sec-1.  Units  of  hardness  and  toughness 
used  to  cai.  late  the  abscissae  of  this  figure  were  Pa  and  Pam1'* ,  respectively. 
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model  of  erosion,  dimensional  analysis  provides 
an  operative  equation  to  describe  erosion  in  terms 
of  dimensionless  groups  of  variables  and  empirical 
constants  determined  from  a  regression  analysis  of 
experimental  data.  The  empirical  constants  can  be 
compared  with  those  determined  from  the  two 
erosion  models  to  yield  an  unbiased  comparison  of 
theory  and  experiment. 

To  apply  dimensional  analysis,  we  assume  that 
the  parameters  that  control  erosion  are  those  given 
in  Equations  4  and  S.  The  volume  loss  per  particle 
impact,  W,  is  then  a  function  of  these  parameters: 

W  =  F(v0,R,p,Ke,H)  (6) 

Using  the  standard  methods  of  dimensional 
analysis  [39],  the  following  functional  relation 
is  obtained: 

W/R3  m  A(K7/RH7f(pv7/H)b  (7) 

where  A  is  a  dimensionless  constant*. 

Because  Young’s  modulus  has  been  suggested 
as  a  variable  that  contributes  to  crack  formation 
during  hardness  indentations  [40],  the  treatment 
given  above  was  extended  to  include  Young’s 
modulus,  E.  The  equation  obtained  is  similar  to 
Equation  7,  but  contains  an  extra  dimensionless 
term,  E/H 

W/R3  =  A'  (Kl/RH'f  (p  t ?/H)b  ( E/Hf 

(8) 

Equation  7  contains  three  dimensionless  groups 
each  of  which  has  physical  meaning.  The  first 
represents  the  ratio  of  the  volume  loss  during 
impact  to  the  volume  of  the  impacting  particle. 
All  other  parameters  being  constant,  the  erosion 
rate  will  increase  as  the  particle  volume  (i.e. 
particle  size)  increases.  The  second  group 
( Kl/RH *)  can  be  represented  as  the  ratio  of  the 
inverse  of  target  brittleness  to  the  size  of  the 
impacting  particle.  This  interpretation  follows 
from  the  fact  that  ( KJH )7  is  a  measure  of  the 
relative  resistance  of  a  target  to  fracture  during  an 
impact  event:  the  higher  the  value  of  (KJH)7, 
the  more  resistant  the  target  will  be  to  fracture. 
The  parameter  (KJH)7  can  be  thought  of  as 
representing  a  critical  scaling  dimension  above 
which  fracture  occurs  during  contact.  The  inverse 
quantity ,  (H/Ke)7 ,  is  a  useful  index  of  “brittleness”. 


[40]  The  third  group,  (pi?/H)  represents  the  ratio 
of  the  particle  energy  density,  i.e.  kinetic  energy 
per  particle  volume,  pi? ,  to  the  hardness,  which 
can  be  considered  as  a  deformation  energy  density. 
The  fourth  dimensionless  constant  in  Equation  8 
can  be  considered  as  the  ratio  of  the  elastic  to  the 
plastic  energy  density. 

The  constants  a  and  b  for  the  dynamic  model 
of  erosion,  have  values  of  —0.67  and  1.58,  res¬ 
pectively.  For  the  quasi-static  erosion  model,  the 
values  of  a  and  b  are  —  0.67  and  1 .22,  respectively. 
The  value  of  c  in  Equation  8  is  zero  for  both 
models.  Hence,  the  two  models  differ  only  in  the 
exponent  of  the  third  dimensionless  group. 

Empirical  values  of  the  constants,  a,  b  and  c  for 
Equation  7  and  8  were  obtained  by  a  multiple 
regression  analysis  of  the  data  (MgO  excluded) 
reported  in  Fig.  2.  The  results  of  the  analysis  are 
given  in  Table  IV.  The  statistics  in  this  table  give 
useful  information  on  the  relative  importance  of 
the  constants  a,  b  and  c  with  regard  to  the  fit  of 
the  erosion  data.  Virtually  the  same  values  of  a 
and  b  and  their  standard  errors  are  obtained 
regardless  of  whether  two  or  three  independent 
variables  are  used  for  the  regression  analysis. 
The  standard  errors  for  a  and  b  are  relatively 
small  (7  and  11%  of  the  mean,  respectively)  and 
the  values  of  ft  computed  for  these  constants  are 
large  and  hence  significant  for  any  reasonable 
level  of  probability.  By  contrast,  the  standard 
error  for  c  (55%  of  the  mean)  is  large,  and  the 
value  of  t  obtained  for  this  constant  is  not  signifi¬ 
cant  at  the  95%  level,  which  suggests  that  the 
value  of  c  reported  in  Table  IV  does  not  differ 
significantly  from  zero.  From  this  discussion  we 
conclude  from  our  results  that  the  wear  rate  does 
not  depend  in  any  significant  way  on  the  ratio  of 
the  Young’s  modulus  to  the  hardness,  E/H* .  This 
conclusion  is  supported  by  the  fact  that  r 1 ,  which 
gives  the  fraction  of  the  variance  accounted  for 
by  the  regression  analysis,  only  changes  from  94 
to  95%  when  E/H  is  added  as  an  independent 
variable.  Consequently,  the  values  of  a  and  b 
determined  from  the  two  parameter  regression 
analysis  will  be  used  for  purposes  of  further 
discussion  in  this  paper. 

As  can  be  seen  from  Table  IV,  the  empirical 
value  for  a,  —0.932,  is  greater  in  absolute  value 


’This  equation  it  derived  in  Appendix  A. 

^For  a  definition  of  t,  tee  any  standard  statistics  text,  e.g.  [41 ). 

*This  conclusion  must  be  tempered  by  the  fact  that  E/H  only  varied  by  a  factor  of  about  2  in  the  present  study.  A 
larger  variation  of  this  parameter  might  indicate  a  significant  dependence  of  wear  on  E!H. 
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TABLE  IV  Determination  of  the  exponents  of  Equations  7  and  8  by  a  multivariable  regression  analysis  (room  tem¬ 
perature  data) 


Exponent 

Regression 

coefficient 

Standard  error 
of  coefficient 

Computed 

t 

Equation  7 

a 

-0.932 

0.110 

-  8.44 

b 

1.384 

0.093 

14.92 

Intercept,  In  A 

-  11.40 

Multiple  correlation 

0.971 

T 3 

0.942 

Standard  error  of  estimate 

0.415 

Equation  8 

a 

-0.905 

0.106 

-  8.56 

b 

1.312 

0.096 

13.66 

c 

-0.669 

0.364 

-  1.84 

Intercept,  In  A 

-  9.84 

Multiple  correlation 

0.975 

r 3 

0.950 

Standard  error  of  estimate 

0.415 

than  the  theoretical  value  of  a,  —0.667,  given  by 
Equations  4  and  5.  The  value  for  b,  1.38,  lies 
approximately  half-way  between  the  value  of  1.22 
predicted  by  the  quasi-static  theory  and  the  value 
of  1.58  predicted  by  the  dynamic  theory.  Using 
the  values  of  a  and  b  from  the  multiple  regression 
analysis.  Equation  8  can  be  expressed  in  a  form 
that  is  similar  to  that  of  Equations  4  and  5 

W<xv2*R3V*K-e'9H04t  (9) 

The  most  significant  difference  between 
Equations  4  and  5  and  Equation  9  is  the  depen¬ 
dence  of  the  wear  rate  on  the  fracture  toughness 
and  the  hardness.  The  exponent  of  Kc  suggests  a 
stronger  dependence  on  this  value  than  is  predicted 
theoretically.  As  Kc  of  the  target  is  increased,  the 
difficulty  of  removing  material  by  chipping 
increases  more  rapidly  than  predicted  by  either 
theoretical  treatment.  Possible  sources  of  this 
variation  involve  the  effect  of  microstructure  on 
erosion  and  the  random  nature  of  the  impact 
process.  These  sources  of  variation  imply  that  the 
models  suggested  to  explain  erosion  may  be  too 
simple  to  account  fully  for  the  effect  of  fracture 
toughness  on  the  erosion  rate .  The  effect  is  micro- 
structure  and  the  random  nature  of  the  particle 
impact  process  will  be  discussed  more  fully  in  a 
later  section  of  this  paper. 

In  view  of  the  fact  that  most  theories  of  erosion 
predict  a  decrease  in  the  erosion  rate  as  the  hard¬ 
ness  in  increased,  the  positive  exponent  of  the 
hardness  in  Equation  9  requires  some  rationaliz¬ 


ation.  An  explanation  for  the  positive  exponent 
for  hardness  in  Equation  9  can  be  developed  from 
a  clow:  examination  of  the  quasi-static  theory  of 
erosion.  In  this  theory,  hardness  determines  both 
the  depth  of  penetration  and  the  maximum  load 
during  impact.  In  the  expression  for  maximum 
load,  hardness  enters  the  equations  with  a  positive 
exponent,  such  that  for  a  fixed  impact  energy  the 
maximum  impact  load  increased  as  the  hardness  is 
increased.  Since  the  amount  of  chipping  is  pro¬ 
portional  to  the  maximum  load  during  impact, 
the  relation  between  load  and  hardness  suggests 
that  the  erosion  rate  increases  as  the  hardness 
increases.  Penetration  is  also  important  because 
it  determines  the  depth  beneath  the  surface 
where  lateral  cracks  form:  the  deeper  the  pen¬ 
etration  (lower  hardness)  the  greater  the  erosion 
rate.  In  the  final  erosion  equation,  the  pentration 
term  and  the  load  term  oppose  one  another  with 
regard  to  hardness  and  the  larger  of  the  two 
determines  the  exponent  for  the  hardness.  In 
Equation  5  the  load  term  dominates,  and  the 
exponent  for  the  hardness  is  positive.  Using  this 
same  line  of  reasoning,  the  results  of  the  regression 
analysis  suggest  that  hardness  effects  surface 
load  to  a  greater  extent  than  penetration  depth, 
resulting  in  a  positive  exponent  for  hardness  in 
Equation  9.  The  fact  that  the  exponent  in  Equation 
9  is  greater  than  0.1 1  suggests  that  the  effect  of 
the  surface  load  term  on  erosion  is  greater  than 
that  predicted  by  the  quasi-static  theory  of 
erosion. 
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4.4.  Microstructural  analysis  of  impact 
damage 

As  noted  in  the  previous  section  of  this  paper, 
the  large  exponent  of  Kc  in  Equation  9  may  have 
its  origin  in  effects  due  either  to  target  micro- 
structure  or  to  the  random  nature  of  the  particle 
impact  process.  The  two  theories  of  erosion 
discussed  in  this  paper  are  predicated  on  the 
assumption  that  the  particles  impact  on  sharp 
comers  and  that  the  type  of  damage  is  similar, 
regardless  of  the  properties  of  the  target  material. 
If  either  assumption  is  not  valid,  then  the  depen¬ 
dence  of  erosion  on  Ke  will  differ  from  that  given 
by  Equations  4  and  S. 

Examination  of  surfaces  that  have  been  impacted 
by  small  numbers  of  particles  yields  information 
on  both  the  type  of  damage  that  occurs  during 
erosion  and  the  relative  number  of  particles  that 
result  in  chipping  from  the  target  surface.  When 
particles  impact  the  target,  they  either  leave 
shallow,  plastic  impressions,  or  small  chipped 
regions  at  the  point  of  impact  (Fig.  6).  The  plastic 
impressions  are  probably  left  by  particles  that 
were  not  oriented  to  impact  on  a  sharp  corner,  but 
on  a  side,  or  edge.  The  residual  plastic  impressions 
in  the  target  surface  suggest  that  the  deformation 
was  not  concentrated  sufficiently  to  nucleate  and 


propagate  surface  cracks,  i.e.  stresses  at  the  impact 
site  did  not  exceed  the  threshold  for  fracture  at 
these  shallow  impressions.  Consequently,  only  a 
fraction  of  particles  that  impact  the  target  surface 
are  effective  in  the  removal  of  material.  If  this 
fraction  depends  on  Kc ,  then  the  erosion  of  the 
target  will  also  depend  on  Kc>  but  in  a  way  not 
given  by  Equations  4  and  5. 

In  our  study  of  the  morphology  of  eroded 
surfaces  we  have  observed  that  as  the  toughness 
of  the  material  is  increased,  the  relative  number 
of  impacts  that  result  in  chipping  is  reduced, 
regardless  of  impact  velocity.  This  observation 
of  impact  behaviour  is  illustrated  in  Fig.  6,  where 
erosion  surfaces  of  glass,  sapphire,  and  silicon 
nitride  are  compared.  As  can  be  seen,  the  fraction 
of  impacts  that  results  in  fracture  and  material 
removal  increases  as  the  fracture  toughness  of  the 
target  decreases.  For  glass,  every  impact  site  in 
Fig.  6a  has  resulted  in  crack  formation.  By  con¬ 
trast,  both  the  sapphire  and  the  silicon  nitride 
have  several  impact  sites  where  plastic  impressions 
were  left,  but  where  crack  formation  was  not 
apparent.  This  observation  suggests  that  the 
functional  dependence  of  erosion  rate  on  Kc 
for  brittle  materials  is  not  completely  described 
by  either  theoretical  treatment  of  erosion 
(Equations  4  and  S),  but  instead  depends  on 


factors  that  are  related  to  the  nucleation  of 
cracks  at  the  impact  site.  Apparently  crack 
nucleation  is  relatively  easier  when  a  “blunt” 
impact  occurs  in  glass  or  silicon,  than  when  it 
occurs  in  the  hot -pressed  materials  used  in  the 
present  study.  Hence,  the  rather  large  dependence 


Figure  6  Single  particle  impact  damage  in  ceramic 
materials.  Optical  micrographs  of  (a)  soda-lime- silica 
glass,  90 m  sec'1 ;  (b)  sapphire,  90 m  sec" 1 ;  (c)  hot-pressed 
silicon  nitride,  90  m  sec" 1 . 


of  erosion  rate  on  Kc  and  H  leported  in  this  paper 
can  be  attributed  in  part  to  statistical  effects  of 
particle  orientation  during  erosion  and  the  ease 
with  which  cracks  nucleate  in  the  target  surface. 

A  second  possible  explanation  for  the  obser¬ 
vation  of  a  larger  than  expected  dependence  of 
erosion  rate  on  Kc  and  H  has  to  do  with  the 
geometry  of  the  cracks  that  form  during  the 
erosion  process.  The  theories  that  have  been 
proposed  to  explain  erosion  assume  that  cracks 
propagate  from  the  impact  site  in  a  self-similar 
fashion,  i.e.  the  cracks  formed  during  impact  are 
geometrically  similar.  Once  nucleated,  the  cracks 
are  assumed  to  propagate  to  the  target  surface. 
Thus,  material  is  removed  from  the  target  by  each 
impact.  In  contrast  to  these  expectations,  micro¬ 
scopic  examination  of  the  target  surface  indicates 
that  the  effectiveness  of  material  removal  from  the 
target  seems  to  depend  on  the  fracture  toughness 
of  the  target.  Thus  for  the  hot-pressed  materials 
used  in  this  study,  cracks  are  often  observed  to 
arrest  within  the  solid  Fig.  7b  and  c  so  that  crack¬ 
ing  during  impact  does  not  result  in  material  loss. 
A  second  or  third  impact  in  the  vicinity  of  the 


primary  impact  site  is  needed  for  material  to  be 
removed  front  the  target.  In  contrast  to  this 
behaviour,  complete  chipping  from  the  primary 
impact  site  is  a  more  frequent  occurrence  for  the 
more  brittle  materials  such  as  silicon  or  glass 
(Fig.  7a).  Thus,  as  the  fracture  toughness  of  the 
material  increases,  the  efficiency  of  material 
removal  per  impact  event  is  less  than  predicted 
theoretically,  and  the  effect  of  Kc  on  the  erosion 
rate  is  greater  than  predicted  theoretically. 

Before  turning  from  the  subject  of  micro¬ 
structure,  it  is  worth  commenting  on  the  erosion 
results  obtained  for  aluminium  oxide  and  sapphire. 
As  can  be  seen  from  Fig.  2,  the  erosion  rates  of 
sapphire  and  sintered  aluminium  oxide  are  approxi¬ 
mately  three  times  that  of  the  hot-pressed  alu¬ 
minium  oxide.  The  difference  in  behaviour  of  the 
two  polycrystalline  materials  is  attributable  to  the 
difference  in  grain  size  of  the  two  materials.  The 
grain  size,  3  to  4pm,  of  the  hot-pressed  material 
was  considerably  smaller  than  the  size  of  the  lateral 
cracks  that  were  formed  upon  impact.  As  a  con¬ 
sequence,  lateral  cracks  interact  with  many  grains 
during  propagation,  and  the  effective  value  of 
Ke  resisting  the  growth  of  lateral  cracks  is  that 
typical  of  polycrystalline  aluminium  oxide, 
~4MPam1/J.  By  contrast,  lateral  cracks  formed 
in  the  sintered  aluminium  oxide,  grain  size  ~  30 
pm,  are  usually  contained  within  a  single  grain, 
and  the  effective  Kc  resisting  crack  growth  is 
more  typical  of  values  obtained  from  single 
crystal  fracture  measurements,  ~  2  MPamI/2 . 


Figure  7  Single  particle  impact  damage  in  ceramic 
materials.  Scanning  electron  micrographs  of  (a)  soda- 
lime-silica  glass,  54  msec'1;  (b)  hot-pressed  aluminium 
oxide,  90m  sec'1;  (c)  hot-pressed  silicon  nitride,  90  m 
sec'1. 
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Using  Equation  9,  the  erosion  rate  for  the  large 
grain  aluminium  oxide  should  be  approximately 
3.7  times  that  obtained  for  the  fine  grained 
aluminium  oxide.  In  Fig.  2,  the  erosion  rate  for 
the  sintered  aluminium  oxide  is  approximately 
33  times  that  of  the  hot -pressed  material,  which 
is  close  to  the  expected  value.  The  fact  that  the 
erosion  rate  of  the  sintered  aluminium  oxide 
is  close  to  that  obtained  for  sapphire  lends  further 
support  for  this  interpretation  of  the  data. 

4.5.  Erosion  at  elevated  temperatures 
As  can  be  seen  by  comparing  Figs.  3  and  4  with 
Fig.  2,  the  temperatures  employed  in  the  present 
study  have  a  marginal  effect  on  the  rate  of  erosion. 
This  finding  is  consistent  with  that  reported  earlier 
by  the  authors  from  a  more  limited  set  of  data 
collected  on  some  of  the  same  materials  studied 
in  this  paper  [7],  Since  dislocation  mobility  is 
enhanced  by  increasing  the  temperature,  it  was 
expected  that  both  the  hardness  and  fracture 
toughness,  and  hence  the  erosion  rate,  would 
be  modified  by  increasing  the  temperature.  Indeed, 
when  loads  are  applied  slowly,  both  the  hardness 
and  toughness  of  ceramic  materials  are  strdngiy 
dependent  on  temperature  [19,  26,  29,  42,  43J. 
The  fact  that  significant  changes  in  the  erosion 
behaviour  are  not  observed  at  elevated  temper¬ 
atures  suggests  that  for  conditions  of  dynamic 
loading,  both  the  hardness  and  the  toughness  are 
invariant  with  temperature.  This  supposition  is 
supported  by  dynamic  toughness  measurements 
on  hot-pressed  silicon  nitride  [44,  45],  and  by  the 
fact  that  cracks  are  observed  to  form  in  soda- 
lime-silica  glass  at  temperatures  above  the  soften¬ 
ing  point  of  this  glass  [46, 47). 

Although  temperature  does  not  play  a  dominant 
role  in  the  erosion  of  ceramics  under  the  conditions 
used  in  the  paper,  minor  differences  between  low 
and  elevated  temperature  behaviour  can  be  attri¬ 
buted  to  plastic  flow.  Because  lateral  cracks  form 
after  the  impact  event  and  are  driven  by  residual 
stresses  at  the  impact  site,  relaxation  of  those 
stresses  or  modification  of  the  resistance  of  the 
target  to  crack  growth  as  a  result  of  plastic  defor¬ 
mation  can  alter  the  size  of  chips  that  are  formed 
after  impact.  Such  effects  are  feasible  when  the 
relaxation  time  of  the  material  for  plastic  flow 
is  less  than  the  time  required  for  the  lateral  cracks 
to  complete  their  growth.  In  an  earlier  study  on 
soda-lime-silica  glass  at  500°  C,  the  temperature 
dependence  of  the  erosion  rate  was  attributed  to 


such  plastic  relaxation  [46],  In  the  present  study, 
the  small  differences  between  low  and  elevated 
temperature  behaviour  may  also  be  attributable 
to  the  same  types  of  processes. 
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Appendix  A: 

Dimensionless  erosion  equation 
This  Appendix  is  written  for  those  readers  who  are 
unfamiliar  with  the  technique  of  dimensional 
analysis.  A  concise  description  of  the  technique 
can  be  found  from  Kay  and  Nedderman  [48]. 
The  procedure  outlined  in  this  reference  is  followed 
here. 

Starting  with  Equation  6,  we  assume  that  the 
wear  rate  can  be  expressed  as  a  power  series 
expansion  of  the  parameters  v0,R,p,Ke  and  H 

*  =  I  **<,#•<)  (Al) 

i 

otj  being  a  dimensionless  coefficient  for  each  term, 
of  the  senes.  The  dimensions  of  each  term  in  the 
expansion  must  equal  the  dimensions  of  W  in  order 
for  Equation  Al  to  be  dimensionally  consistent. 

When  the  dimensions  are  substituted  for  the 
parameters  in  Equation  Al,  each  term  in  the 
expansion  must  have  the  following  dimensional 
form: 

L3  =  (L/T)“tL6'(M/L3)c<(M/T2L,/J)d'(M/T1L)*' 

(A2) 

where  L,  T,  and  M  represent  the  dimensions  of 
length,  time,  and  mass. 

Equating  the  exponents  for  each  dimension 
we  obtain  three  simultaneous  equations  in  terms 
of  the  exponents  a,,  b,,  ct,  dt,  and  e,.  If  two  of  the 
unknowns  are  selected  as  independent  variables, 
the  other  unknowns  can  be  expressed  in  terms 
of  these  two  variables.  For  example,  if  ct  and  d, 
are  selected  then  the  following  equations  are 
obtained  for  at,  bt,  and  et: 

a,  =  —  3  —  2c( 

b ,  =  3 -d,l 2 

e,  =  ~c,-d,  (A3) 

If  these  are  substituted  into  Equation  Al.  the 
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following  expression  is  obtained  for  the  wear  rate 
W: 

W/R 3  =  £  a, 

(A4) 

Since  the  two  theories  developed  to  explain 
erosion  are  power  functions  of  the  variables 
given  in  Equation  A4,  only  one  term  in  the  series 
need  be  retained  in  order  to  compare  the  dimen¬ 
sional  analysis  with  the  theoretical  expressions 
given  by  Equations  4  and  5.  Hence,  the  following 
relation  is  obtained  for  the  erosion  rate: 

W/R3  =  A  {Kl/RH2?  [pVl/H)b  (AS) 

where  A  is  a  dimensionless  constant  and  the 
exponents  of  Equation  A4  have  been  written  as 
a  and  b.  Equation  AS  is  identical  to  Equation  7 
of  the  text.  As  noted  earlier,  the  undetermined 
constants  a  and  b  are  evaluated  by  an  empirical 
fit  of  erosion  data. 

Equation  A5  is  not  a  unique  dimensionless 
representation  of  the  parameters  that  control 
erosion.  For  example,  if  bt  and  c,  had  been  selected 
as  the  independent  constants,  then  the  following 
erosion  equation  would  have  been  obtained 

W(H/KC)6  =  A  '(v2p/H)c  (RH2/K3)b  (A6) 


Equations  A5  and  A6  can  be  shown  to  be  equiv¬ 
alent  by  dividing  both  sides  of  the  equation  by 
(RH2/K%)3.  By  systematically  solving  for  all 
possible  combinations  of  the  exponents  in  Equation 
Al,  five  variants  of  Equation  A5  were  found. 
These  could  all  be  reduced  to  Equation  AS  by 
judicial  manipulation  (multiplying  or  dividing) 
of  the  dimensionless  parameters  Kl/RH3  and 
pit/H. 

There  is  a  certain  arbitrariness  in  selecting 
one  of  the  dimensionless  equations  for  a  com¬ 
parison  with  the  experimental  data.  We  justify 
the  selection  of  Equation  A5  on  the  basis  of  its 
simple  form  and  the  ease  with  which  the  dimen¬ 
sionless  variables  K3/RH 2  and  ptf/H  can  be  given 
physical  interpretation.  Furthermore,  this  arrange¬ 
ment  of  the  variables  in  Equation  A5  separates 
the  variants  used  in  the  present  study  more  effec¬ 
tively  than  the  others,  and  permits  us  to  compare 
the  theoretical  equations  with  the  results  of  the 
dimensional  analysis  more  readily.  Regardless 
of  which  form  of  the  dimensionless  analysis  is 
used,  one  can  show  that  they  are  all  equivalent, 
provided  the  error  is  minimized  in  the  term 
containing  the  wear  rate,  i.e.  W/R3.  This  equiv¬ 
alence  can  be  demonstrated  by  using  the  basic 
equations  for  a  multiple  regression  analysis  [41]. 
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T  ABLE  IB  Summary  of  erosion  data 

Material 

Temperature  (°  C) 

Particle  velocity 
(m  sec'1 ) 

Erosion  rate 
(mm3) 

Hot-pressed  Si3N4 

25 

94 

7.4 

X 

IO'7 

(0.8)* 

73 

4.0 

X 

IO'7 

(0.6) 

37 

9.9 

X 

10-" 

(2.6) 

500 

125 

1.8 

X 

10‘‘ 

(0.2) 

90 

8.1 

X 

10-’ 

(2.0) 

54 

2.2 

X 

IO'7 

(0.4) 

1000 

125 

2.5 

X 

10-‘ 

(0.2) 

90 

9.9 

X 

IO'7 

(0.5) 

Hot-pressed  SiC 

2S 

94 

7.2 

X 

10- ' 

(1.5) 

73 

5.3 

X 

10-7 

(0.2) 

37 

1.4 

X 

10-7 

(0.3) 

Hot-pressed  Al,03 

25 

94 

1.1 

X 

10'‘ 

(0.2) 

73 

6.0 

X 

10-7 

(0.6) 

37 

1.3 

X 

io-7 

(0.4) 

500 

125 

2.5 

X 

10  ‘ 

(0.1) 

100 

1.5 

X 

io-* 

(0.1) 

54 

4.1 

X 

10' 

(0.1) 

1000 

125 

2.5 

X 

IO'* 

(0.1) 

81 

8.7 

X 

10'7 

(0.6) 

54 

3.8 

X 

IO'7 

(0.4) 
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TABLED  continued 


Material 


Sintered  At,  O 


Temperature  (*  Q 


Particle  velocity 
(m  tec"1 ) 


Sintered  A1,0 


Sintered  MgO 


Fused  Silica 


Soda-lime - 
silica  glass 


Erosion  ra 
(mm*) 


3.2 

X 

10- 

(0.4) 

2.0 

X 

10- 

(0.1) 

4.3 

X 

10- 

(1.3) 

9.1 

X 

10* 

(1.0) 

3.9 

X 

10- 

(1.1) 

8.8 

X 

10- 

(3.6) 

6.7 

X 

10- 

(2.1) 

2.8 

X 

10- 

(0.7) 

1.8 

X 

10- 

(0.6) 

7.0 

X 

10- 

(2.5) 

4.6 

X 

10- 

(0.4) 

1.2 

X 

10- 

(0.1) 

5.5 

X 

10- 

(0.1) 

9-1 

X 

10- 

(1.0) 

3.9 

X 

10- 

(1.1) 

8.8 

X 

10- 

(3.6) 

9.8 

X 

10- 

(0.7) 

3.8 

X 

10- 

(0.3) 

6.4 

X 

10- 

(0.6) 

2.1 

X 

10- 

(0.1) 

4.1 

X 

10- 

(0.3) 

1.4 

X 

10- 

(0.1) 

6.5 

X 

10- 

(0.7) 

2.6 

X 

10- 

(0.2) 

7.7 

X 

10- 

(1.0) 

4.6 

X 

10- 

(0.7) 

7.3 

X 

10- 

(0.7) 

9.8 

X 

10- 

(0.2) 

2.3 

X 

10- 

(0.3) 

1.6 

X 

10- 

(0.2) 

6.0 

X 

10- 

(0.5) 

4.7 

X 

10- 

(0.0) 

3.2 

X 

10- 

(0.7) 

1.7 

X 

10- 

(0.4) 

7.1 

X 

10- 

(0.7) 

3.2 

X 

10- 

(0.7) 

4.4 

X 

10- 

(1.2) 

2.2 

X 

10- 

(0.6) 

8.3 

X 

10- 

(0.4) 

2.3 

X 

10- 

(0.5) 

The  numbers  in  parentheses  are  the  standard  deviation  of  each  fit.  Two  to  ten  erosion  measurements  were  used  to 
determine  each  erosion  rate. 
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B.R.  Lawn  and  S.M.  Wiederhorn 


Fracture  and  Deformation  Division 
National  Bureau  of  Standards 
Washington,  D.C.,  U.S.A. 


ABSTRACT:  The  nature  of  contact- induced  surface  damage  in  brittle  materials, 
and  the  fracture  mechanics  principles  used  to  describe  this  damage, 
are  surveyed.  The  importance  of  understanding  the  elastic  and 
plastic  deformation  processes  which  precede  fracture  is  emphasized. 
Strength  and  erosive  wear  properties  are  intimately  connected  to 
the  contact  damage  mechanics. 

RESUME:  La  nature  de  la  deterioration  par  contact  de  la  surface  des 

materiaux  fragiles  et  les  principes  de  mecanique  des  fractures 
utilises  pour  decrire  cette  deterioration  sone  rappeles.  On 
souligne  1' importance  de  la  comprehension  des  processus  de 
deformation  elastique  et  plastique  qui  prSc&dent  la  fracture. 

La  force  et  les  proprietes  d'erosion  sont  intement  lies  aux 
mecanismes  de  doomage  par  contact. 


INTRODUCTION 


The  ubiquitous  surface  damage  that  characterises  highly  brittle  materials,  notably  glasses 
and  ceramics  (i.e.,  solids  with  covalent/ionic  bonding),  is  due  to  local  stress  concentra¬ 
tions  that  occur  whenever  contact  is  made  with  a  small,  hard  object.  Microfracturc  centers 
that  seriously  degrade  the  strength  are  often  introduced  by  the  processes  used  to  finish 
the  surfaces  (e.g.,  machining) ,  or  by  particle  impingement  incurred  in  subsequent  handling 
and  storage.  Unless  extreme  precautions  are  taken  to  avoid  all  spurious  contact  events 
(e.g. ,  as  is  done  with  the  coating  of  freshly  drawn  optical  fibres  in  dust-free  atmospheres) 
such  degradation  is  generally  inevitable.  Contact  damage  also  holds  the  key  to  the  ero¬ 
sive  wear  and  abrasion  properties  of  brittle  materials.  A  proper  understanding  of  the 
underlying  mechanisms  of  deformation  and  fracture,  using  the  controlled  methods  of  "inden¬ 
tation  fracture  mechanics",  has  accordingly  become  a  major  research  goal  in  the  area  of 
brittle  design. 

In  the  present  paper  the  current  state  of  this  understanding  is  summarized. 

Our  goal  is  not  a  comprehensive  survey  of  the  field:  rather,  we  seek  to  draw  attention 
to  certain  broad  features  of  the  brittle  indentation  problem  that  might  be  considered  to 
bear,  however  indirectly,  on  the  theme  of  this  meeting.  Reference  is  made  to  several 
review  articles  [1-4]  for  those  who  wish  to  pursue  the  subject  in  greater  detail. 
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The  outline  of  our  presentation  is  as  follows.  First  we  define  what  we  mean 
by  a  brittle  solid.  We  then  argue  that  indentation  events  can  be  classified  into  two 
main  types,  "blunt"  or  "sharp",  according  to  whether  the  material  response  to  fracture 
is  essentially  elastic  or  plastic.  In  normal  loading  the  cracks  are  shown  to  have  well- 
defined,  penny-like  geometries;  superposition  of  a  tangential  loading  component  modifies 
these  geometries  significantly.  Fracture  mechanics  relations  are  given  for  some  of  the 
more  important  of  the  crack  geometries.  Finally,  the  role  of  indentation  fracture  des¬ 
criptions  in  formulating  theories  of  strength  and  wear  is  discussed. 

BRITTLE  MATERIALS  IN  RELATION  TO  THE  CONTACT  PROBLEM 

Ideally  brittle  materials  are,  by  definition,  essentially  characterized  by  a  completely 
elastic  response  up  to  the  point  of  fracture.  In  certain  instances  stresses  and  strains 
close  to  the  theoretical  strength  of  the  molecular  structure  can  be  sustained  without 
detectable  signs  of  permanent  deformation.  Coated  silica  glass  fibers,  for  example,  show 
complete  recovery  after  undergoing  tensile  strains  of  up  to  15%.  The  materials  which 
fall  most  readily  into  this  category  are  those  with  large  components  of  covalent  bonding, 
for  which  there  exists  a  strong  intrinsic  resistance  to  shear-activated  deformation 
processes  [5] . 

However,  even  the  most  brittle  of  materials  can,  if  subjected  to  sufficiently 
large  constraining  hydrostatic  compressions  to  inhibit  the  onset  of  fracture,  be  deformed 
irreversibly  [6].  (This  statement  is,  of  course,  a  self-evident  truth  to  those  concerned 
with  the  geomechanical  behaviour  of  rocks.)  Hardness  indentations  provide  us  with  the 
simplest  means  of  demonstrating  this  phenomenon;  the  stress  field  immediately  beneath 
the  contact  area  is  intensely  compressive,  with  a  substantial  component  of  superposed 
shear  [7].  Thus  residual  impressions  can  be  made  on  the  surface  of  any  material,  includ¬ 
ing  diamond,  with  a  suitably  penetrative  indenter. 

Whereas  in  metals  the  nature  of  the  deformation  processes  which  operate  within 
the  contact  zone  is  reasonably  well  understood,  in  some  of  the  more  brittle  materials  the 
analogous  processes  remain  obscure.  It  is  clear  from  the  magnitude  of  the  contact  pres¬ 
sures  that  the  harder,  covalent  structures  are  being  stressed  to  their  theoretical  limit. 
At  this  level  the  classical  descriptions  of  slip  by  dislocation  motion  no  longer  strictly 
apply;  instead  it  becomes  more  useful  to  consider  the  deformation  modes  in  terms  of  an 
extended,  cooperative  breakdown  of  the  structure.  This  is  not  to  say  that  structures 
which  undergo  this  kind  of  deformation  are  incapable  of  being  "dislocated"  by  shear 
processes.  Indeed,  structural  dislocations  have  been  clearly  identified  in  transmission 
electron  microscopy  observations  by  Hockey  at  indentation  sites  in  a  number  of  hard 
crystalline  materials  [8-10].  However,  the  configurations  observed  do  not  always  corres¬ 
pond  to  normal  crystallographic  slip  planes  or  directions.  In  fact,  the  recent  identifi¬ 
cation  of  analogous  shear  processes  in  soda-lime  glass  [11}  would  appear  to  indicate  that 
crystallographic  considerations  are  no  longer  of  primary  importance  in  the  constrained 
deformation  of  this  class  of  solid. 
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A  characteristic  feature  of  the  contact  deformation  zone  in  highly  brittle 
materials  is  its  strong  confinement  to  the  region  immediately  below  the  surface  impres¬ 
sion.  There  is  no  mechanism  for  relaxation  of  the  "plastic"  strains  as  there  is  in  most 
metals,  where  extensive,  long-range  slip  or  twinning  can  usually  occur  without  obstruc¬ 
tion.  Instead,  these  strains  have  to  be  accommodated  elastically  by  the  surrounding 
matrix.  Consequently,  high- intensity  residual  stress  fields  can  develop,  and  these 
fields  can  exert  a  strong  influence  on  subsequent  mechanical  response  of  the  material. 

A  second  characteristic  feature  of  the  contact  process  in  brittle  materials 
is  the  great  ease  with  which  microcracks  initiate  and  propagate.  In  a  covalent  material 
like  silocon,  for  instance,  it  is  almost  impossible  to  produce  crack-free  impressions, 
with  even  the  most  delicate  of  routine  hardness  testing  machines.  In  this  context  it 
may  be  noted  that  a  tensile  stress  component,  however  small  in  comparison  to  the  hydro¬ 
static  compression  within  the  deformation  zone,  is  generally  unavoidable  in  the  matrix 
contact  field  [1]. 

BLUNT  VERSUS  SHARP  CONTACT 

The  nature  of  the  stress  field  in  an  indentation  experiment  depends  strongly  on  the  geo¬ 
metry  of  the  contacting  surfaces,  as  well  as  on  the  mechanical  properties  of  the  materials 
involved.  We  shall  be  working  on  the  premise  here  that  the  indenter  material  is  suffi¬ 
ciently  "hard"  relative  to  that  of  the  test  piece  to  be  effectively  rigid.  It  is  then 
convenient  to  distinguish  two  extreme  types  of  indentation  field  [1] :  "blunt",  in  which 
the  contact  pressure  increases  monotonically  with  load  such  that  the  deformation  prior 
to  fracture  is  completely  elastic;  "sharp",  in  which  the  contact  pressure  is  in  excess  of 
that  required  to  produce  irreversible  deformation  at  all  stages  of  loading.  The  intensity 
of  the  stress  field  in  the  former  case  is  controlled  by  the  elastic  moduli,  in  the  latter 
by  a  combination  of  the  elastic  moduli  and  the  hardness. 

Blunt  Contact 

The  classical  example  of  the  blunt  contact  is  the  Hertzian  stress  problem,  discussed  at 
some  length  by  others  in  this  volume.  Experimentally,  the  Hertzian  stress  field  is  most 
simply  generated  by  pressing  a  sphere  onto  a  semi-infinite  solid.  Solutions  for  the 
field  are  obtained  from  the  basic  equations  of  linear  elasticity.  Essentially,  the 
normal  stresses  are  all  highly  compressive  in  a  drop-shaped  zone  immediately  below  the 
indenter,  but  become  moderately  tensile  at  and  outside  the  contact  circle.  These  tensile 
stresses  are  extremely  inhomogeneous  in  this  near-contact  region,  falling  off  dramatically 
along  subsurface  stress  trajectories  (1,12].  Remote  from  the  contact  zone  the  stress 
field  tends  asymptotically  to  the  corresponding  Boussinesq  field  for  concentrated  point 
loading  [1,13]. 

At  a  critical  load  in  the  Hertzian  contact  a  well-defined,  cone-shaped  crack 
"pops  in"  from  the  specimen  surface.  Figure  1  shows  such  a  Hertzian  fracture  in  glass. 


The  stresses  at  the  surface  trace  prior  to  pop-in  are  generally  well  below  the  theoretical 
limiting  strength  of  the  material  structure,  indicating  that  initiation  must  occur  from 
pre-present  flaws.  The  crack  first  runs  from  the  critical  surface  flaw  into  a  shallow 
ring  just  outside  the  contact  circle,  then  propagates  downward  into  its  characteristic 
cone  geometry  until  sufficiently  remote  from  the  loading  center,  at  which  point  it  becomes 
highly  stable. 


Figure  1  -  Cone  Crook  in  Glaee ;  Ba.ee  Diameter  30  rrm,  Indentation 
Load  40  MN.  After  [J5J. 


The  mechanics  of  formation  of  Hertzian  cone  cracks  is  complicated  by  the 
extreme  inhomogeneity  of  the  near-contact  stress  field  through  which  the  growth  occurs. 

To  ignore  this  inhomogeneity  and  assume  that  instability  ensues  when  the  surface  stresses 
reach  the  tensile  strength  of  the  material  is  to  overlook  the  essence  of  the  general 
contact  fracture  phenomenon.  In  accordance  with  modern-day  fracture  mechanics  procedure 
it  is  necessary  to  compute  a  "stress  intensity  factor",  representing  the  driving  force 
for  the  fracture,  as  an  integral  of  actual  stresses  (weighted  with  an  appropriate  Green's 
function)  over  the  prospective  crack  path.  The  first  such  analysis  was  carried  out  by 
Frank  and  Lawn  [12],  who  showed  that  the  critical  load  for  cone-crack  pop-in,  under  equili 
brium  conditions  of  fracture,  is  given  by 

Pc  -  A,K2r/E  ,  (1) 


where  r  is  the  sphere  radius,  K,  is  the  critical  stress  intensity  factor  for  crack  exten¬ 
sion  (i.e.,  the  material  "toughness",  characterizing  the  intrinsic  resistance  to  fracture) 
E  is  Young's  modulus  and  A*  is  a  dimensionless  constant.  This  result  provided  the  first 
analytical  derivation  of  the  long-standing  empirical  "law"  of  Auerbach,  1891,  Pc  «  r  [14]. 
Considerable  interest  has  been  shown  in  this  law  because  the  Hertzian  equations  for  the 
surface  stresses  give  o  *  P/r2  so  that,  in  combination  with  equation  (1),  we  obtain 
ac  «  1/r;  i.e.,  the  simplistic  concept  of  a  critical  stress  criterion  for  fracture  is 
clearly  in  violation.  Note  that  this  last  expression  implies  the  suppression  of  crack 
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initiation  at  small  sphere  radii  -  the  possibility  ot'  inducing  precursor  plasticity  thus 
becomes  stronger  for  "sharper"  indenters. 

The  size  to  which  the  cone  crack  grows,  once  initiated,  is  determined  by  the 
far-field  conditions.  Roesler  [IS],  using  a  dimensional  analysis  argument,  showed  that 
the  crack  size  increases  with  load  according  to  c  “  P^3.  A  more  detailed,  fracture 
mechanics  treatment  [16]  gives,  at  equilibrium, 

P/c*2  »  A2Kc  ,  IP  >  Pc)  ,  (2) 

where  A2  is  another  dimensionless  constant. 

It  is  interesting  to  note  that  neither  of  equations  (1)  and  (2)  is  sentitive 
to  the  size  of  the  critical  flaw  from  which  the  cone  crack  initiates.  This  insensitivity 
is  another  characteristic  of  the  general  contact  problem,  although  it  is  implicit  in  the 
derivation  of  the  above  expressions  that  there  is  always  a  sufficiently  high  density  of 
surface  flaws  present  to  guarantee  the  initiation  condition.  It  may  also  be  noted  that 
the  formulations  contain  the  toughness  K^,  the  definition  of  which  strictly  implies  a 
configuration  of  mechanical  equilibrium  (unstable  in  relation  to  equation  (1),  stable  in 
equation  (2)).  In  practice,  brittle  materials  are  susceptible  to  chemically-enhanced 
rate-dependent  crack  growth,  particularly  in  the  presence  of  atmospheric  moisture,  in 
which  case  equation  (1]  tends  to  overestimate  the  critical  load  and  equation  (2)  to 
underestimate  the  crack  size. 

The  mechanics  of  fracture  are  modified  somewhat  by  the  superposition  of  a  tan¬ 
gential  force  component  onto  the  normal  loading  configuration.  Stress  solutions  for  the 
case  of  complete  slippage  at  a  sliding  contact  interface  have  been  described  by  Hamilton 
and  Goodman  [17].  The  main  effect  of  the  tangential  force  is  to  enhance  the  tensile 
stresses  at  the  trailing  edge  of  the  advancing  sphere,  thereby  destroying  the  axial 
symmetry  of  the  field.  Consequently,  the  critical  load  for  fracture  is  reduced,  drama¬ 
tically  at  the  higher  coefficients  of  friction,  and  the  cracks  form  only  partial  (highly 
distorted)  cones  [18].  The  modified  field  for  a  coefficient  of  friction  0.1  is  shown  in 
Figure  2,  together  with  a  micrograph  for  the  corresponding  conditions  in  a  sliding  fric¬ 
tion  test  on  a  glass  surface  [19].  Formulations  analogous  to  equations  (1)  and  (2)  for 
the  sliding  sphere  have  been  attempted  by  several  workers,  but  a  strong  sensitivity  to 
starting  assumptions  in  the  analysis  has  led  to  considerable  divergence  in  the  various 
predictions,  particularly  in  the  critical  load. 

Sharp  Contact 

We  alluded  in  the  previous  subsection  to  the  increasing  prospect  of  precursor  inelastic 
deformation  with  diminishing  indenter  radius.  In  the  limit  of  zero  radius,  e.g.,  as 
with  an  ideally  sharp  Vickers  or  Knoop  pyramidal  indenter,  such  deformation  becomes 
unavoidable;  a  finite  load  cannot  be  supported  by  a  point  contact  without  exceeding  the 
elastic  limit.  The  stress  field  beneath  the  indenter  is  considerably  more  complex  than 
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in  its  Hertzian  elastic  counterpart,  and  simplistic  elastic/plastic  models  have  to  be 
devised  to  provide  a  necessary  framework  for  fracture  mechanics  analysis.  The  simplest 
and  most  widely  adopted  of  the  models  is  that  of  the  expanding  internal  cavity,  in  which 
a  pressurized  spherical  volume  (hardness  impression)  induces  plasticity  in  an  immediate 
annular  surround  volume  (deformation  zone),  the  whole  being  constrained  in  an  infinite 
elastic  matrix  [20,21].  Despite  clear  shortcomings,  e.g.,  failure  to  allow  for  stress 
relaxation  at  the  specimen  free  surface  and  failure  to  match  the  Boussinesq  point-force 
solution  in  the  far  field,  models  of  this  type  are  useful  for  their  amenability  to  closed- 
form  solution. 


Figure  2  -  Sliding  Contact  Fracture  Beneath  Sphere ,  Friction 
Coefficient  0.1.  Schematic  half-surface  and  side 
views  show  contours  of  maximum  tension  (solid  curves) 
and  principal  stress  trajectories  (dashed  curves) 
about  contact  circle  (shaded).  After  [18].  Micrograph 
shows  partial  cone  tract  on  glass ,  width  25  vm,  produced 
by  sphere  of  radius  1.5  mm  at  load  10  0.  After  [19]. 
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The  fracture  pattern  produced  in  sharp  indentation  has  some  geometrical  var¬ 
iants,  but  two  basic  crack  types  are  generally  distinguishable;  median/radial,  and 
lateral,  cracks  [13].  Figure  3  is  an  illustrative  example,  again  in  glass.  It  is  found 
that  the  initiation  of  such  cracks  is  accomplished  with  relative  ease  in  highly  brittle 
materials,  even  on  surfaces  of  the  greatest  perfection,  suggesting  that  the  percursor 
deformation  processes  are  capable  of  producing  their  own  fracture  nuclei  [11].  Each 
median  crack  forms  on  a  symmetry  plane  containing  the  contact  axis  and  an  impression 
diagonal,  and  leaves  a  characteristic  radial  trace  on  the  indented  surface;  as  seen  in 
Figure  3,  more  than  one  crack  of  this  type  can  be  formed  in  a  given  contact,  depending 
on  the  indenter  geometry.  Lateral  cracks  form  in  a  saucer-like  manner,  starting  from 
near  the  base  of  the  deformation  zone  and  spreading  outward  closely  parallel  to  the 
specimen  surface;  in  severe  cases  the  laterals  can  extend  to  the  surface,  causing  chipping 
As  with  the  Hertzian  cones  these  cracks,  once  developed,  are  highly  stable. 


Figure  3  -  Median/ Radial. /Lateral  Fracture  System  for  Vickers  Indenta¬ 
tion  on  Glass;  Radial  Crack  Diameter  1.3  rm ,  load  100  N. 

Courtesy  of  B.J.  Hookey. 

The  conditions  for  the  initiation  of  indentation  cracks  beneath  sharp  indenters 
have  only  recently  been  considered.  Lawn  and  Evans  [22]  used  a  fracture  mechanics  proce¬ 
dure  somewhat  analogous  to  that  in  the  earlier  cone  crack  analysis  [12],  but  with  the 


intensity  of  the  stress  field  now  determined  by  a  load- invariant  hardness  H.  In  their 
scheme  fracture  occurs  when  the  spatial  extent  of  the  field,  as  determined  by  the  con¬ 
tact  dimension,  is  sufficiently  large  to  make  available  flaws  within  the  deformation  tone 
unstable.  Accordingly,  the  critical  load  is  found  to  be 

Pc  -  Aj(E/H)(Kc/H)3Kc  ,  (3) 

where  As  is  a  dimensionless  term  involving  the  characteristic  indenter  half-angle. 

Whereas  the  analysis  is  acknowledged  to  be  deficient  in  its  ability  to  specify  the 
absolute  value  of  Aj,  Equation  (3)  serves  well  as  a  basis  for  ranking  materials.  For 
example,  comparison  of  a  moderately  hard  steel  with  a  typical  silicate  glass  leads  to  a 
predicted  ratio  of  2  103  -  10*  in  Pc  values,  thus  explaining  the  brittleness  of  the 
latter  [7] . 

In  considering  the  mechanics  of  well-developed  median  and  lateral  cracks  we 
can  take  note  of  one  geometrical  feature  common  to  all  far-field  indentation  fracture 
configurations;  their  ultimate  fronts  tend  to  be  circular,  so  they  have  the  essential 
character  of  center- loaded  penny  cracks  [16].  For  the  median  cracks,  which  approximate 
to  half-pennies  centered  on  the  point  of  contact,  detailed  analysis  provides  a  relation 
between  load  P  and  crack  radius  c,  [23] 

P/c3*  -  A^H/E)3*!^  .  (P  >  Pc) 

in  direct  analogy  to  equation  (2),  with  An  again  incorporating  the  indenter  half-angle. 
However,  if  the  final  configuration  of  the  median  crack  bears  a  strong  underlying  resem¬ 
blance  to  that  of  its  cone  crack  counterpart,  its  route  to  this  configuration  during 
the  indentation  cycle  is  of  a  totally  different  nature.  For,  in  addition  to  the  usual 
crack  driving  force  associated  with  the  reversible  (elastic)  component  of  the  stress 
field,  there  is  also  a  contribution  from  the  irreversible  (plastic)  component,  and  it 
is  the  second  of  these  components  which  dominates  in  the  median  crack  evolution.  Thus, 
much  of  the  crack  evolution  occurs  as  the  indenter  is  being  unloaded,  particularly  along 
the  radial  traces  in  the  near-surface  region  [23,24].  Indeed,  in  materials  susceptible 

to  moisture-enhanced  crack  growth,  the  radial  crack  traces  can  be  observed  to  continue 

extending  long  after  the  indenter  has  been  removed. 

Derivation  of  a  similar  expression  for  lateral  cracks  is  even  less  straight¬ 
forward,  owing  to  the  added  complication  of  the  nearby  specimen  free  surface.  Whereas 
for  the  more  penetrative  median  and  cone  cracks  the  free  surface  is  unlikely  to  have  a 
strong  influence  on  the  strain  energy  distributions  which  govern  crack  extension,  the 
same  is  certainly  not  true  of  the  highly  compliant  material  portions  immediately  above 
the  lateral  plane.  Thus,  although  the  lateral  system  retains  the  penny-like  character, 
it  is  more  realistic  to  regard  it  in  terms  of  a  thin,  circular  disc  of  material  built  in 

at  its  circumference  to  a  rigid  matrix  rather  than  in  terms  of  the  usual  embedded  con¬ 

figuration.  Using  the  theory  of  elastic  plates  to  determine  the  energetics  of  this 
system  a  result  analogous  to  equation  (4)  can  be  derived  [25] ,  except  that  now  the 


relation  between  crack  sice  and  load  is  of  the  form  c  *  P^8  and  the  dependence  on  material 
parameters  is  somewhat  more  complex.  Again,  the  bulk  of  the  crack  growth  takes  place 
during  indenter  withdrawal,  re-emphasising  the  crucial  role  of  residual  stresses  when 
inelastic  deformation  accompanies  fracture. 

When  a  sharp  indenter  is  translated  across  the  surface  of  a  brittle  material  it 
leaves  a  scratch  track.  Again,  the  introduction  of  a  frictional  force,  this  time  deter¬ 
mined  by  a  "ploughing”  process,  enhances  the  fracture  development  [26].  In  environment- 
sensitive  materials  rate  variables,  such  as  scratch  velocity,  become  a  factor.  In  accor¬ 
dance  with  the  arguments  leading  to  equation  (3)  it  is  found  that  "light"  loads  tend  to 
produce  smooth  scratches,  "heavy"  loads  fragmented  scratches.  In  the  former  case  the 
action  is  one  of  "polishing",  in  the  latter  of  "abrasion"  [1] .  Figure  4  shows  an  example 
of  a  scratch  made  just  above  the  threshold  load  in  glass.  (Note  that  a  cursory  examina¬ 
tion  of  the  surface  in  this  case  would  create  the  false  impression  of  a  purely  plastic 
contact  event.)  Comparison  of  Figure  4  with  Figure  3  shows  that  the  sliding  motion  tends 
to  suppress  one  set  of  median  cracks,  thus  producing  an  essentially  linear  flaw.  At 
excessively  high  loads  extensive  chipping  occurs  along  the  length  of  the  scratch  as  the 
laterals  turn  upward  to  intersect  the  surface.  This  can  lead  to  physical  removal  of  the 
deformation  zone,  and  thereby  of  the  source  of  the  residual  driving  force  on  the  sub- 
survace  median  crack. 


Figure  4  -  Scratch  on  Glass,  ~  10  \jm  wide,  made  by  sliding  Vickers 
Indenter  at  load  1  N.  After  [26] . 
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PRACTICAL  APPLICATIONS 


Extensive  experimental  confirmation  of  the  preceding  theoretical  formulations  has  pro¬ 
vided  us  with  a  sound  foundation  for  modelling  various  practical  processes  in  brittle 
glasses  and  ceramics.  We  shall  focus  our  attention  on  two  of  the  more  important  of  these 
processes,  strength  degradation  and  erosive  wear.  In  the  first  case  it  is  the  penetrative 
cone  or  median  crack  which  is  the  governing  element;  in  the  second  case  the  lateral 
system  dominates. 

Strength  Degradation 

The  strength  of  a  brittle  material  is  controlled  by  the  size  of  the  largest  flaw  in  the 
surface  or  bulk.  In  many  cases  this  dominant  flaw  results  from  a  well-characterized  con¬ 
tact  event,  e.g.,  from  machining  during  fabrication  or  particle  impact  in  service.  Failure 
then  becomes  a  function  of  the  contact  variables,  which  can  often  be  accurately  specified. 
This  accuracy  can  be  optimized  by  deliberately  introducing  controlled  indentation  cracks 
into  prospective  test  pieces.  Apart  from  circumventing  many  uncertainties  associated 
with  the  usual  statistical  approach  to  the  strength  of  materials  with  unclassified  flaws, 
this  course  opens  up  many  avenues  to  systematic  materials  evaluation  [4] .  Moreover,  by 
observing  how  "macroscopic"  contact- induced  cracks  respond  under  an  applied  tensile 
stress  we  can  gain  valuable  physical  insight  into  the  nature  of  "microscopic"  natural 
flaws  which  are  generally  undetectable. 

As  an  example  of  the  type  of  information  that  can  be  obtained  from  strength 
degradation  studies  we  present  results  from  a  ballistic  impact  stuty  on  glass  surfaces 
using  steel  spheres.  Figure  5  [27].  The  micrograph  shows  the  damage  incurred  at  one 
specified  impact  velocity,  and  the  graph  shows  how  the  remaining  strength  varies  systema¬ 
tically  with  this  velocity.  The  solid  curve  through  the  data  points  is  a  prediction 
from  the  conventional  strength  formula  for  spontaneous  failure  from  flaws  of  characteris¬ 
tic  size  c«. 


S  -  K./YcJ'2 

where  Y  is  a  geometrical  constant  close  to  unity;  eg  is  evaluated  in  accordance  with 
equations  (1)  and  (25  for  blunt  indenters,  using  the  Hertzian  contact  theory  to  eliminate 
load  as  an  unknown  in  favour  of  velocity  [27],  We  note  that  no  dagradation  occurs  in 
Figure  5  below  the  threshold  velocity  for  cone  crack  pop  in,  at  which  point  the  strength 
drops  abruptly;  beyond  the  threshold  the  strength  falls  off  steadily  with  increasing 
velocity. 


Figure  5  -  Strength  Degradation  of  Glass  Impacted  by  Spheres. 

Micrograph  shows  surface  damage  after  impact  with 
35-meah  glass  beads;  diameter  of  cone  c^ack  surface 
trace  400  \jm,  incident  velocity  100  me"  .  Graph 
shows  remaining  strength  as  function  of  impact  velo¬ 
city  with  steel  spheres ,  diameter  0.8  mm.  After  [27]. 

The  strength  behaviour  of  materials  containing  median  flaws  induced  by  sharp 
contact  is  more  complicated,  owing  to  the  active  role  of  residual  stresses  from  the 
plastic  enclave  in  post-contact  crack  growth  [28] .  As  a  result  of  a  strong  stabilizing 
influence  of  the  residual  component  in  the  net  driving  force  for  fracture,  the  medians 
do  not  propagate  spontaneously  to  failure  in  the  manner  of  cone  cracks,  but  rather 
undergo  a  precursor  stage  of  growth  from  Co  to  a  critical  size  cm  before  attaining  an 
instability  configuration.  Some  data  taken  from  direct  observations  of  radial  crack 
growth  during  application  of  a  tensile  stress  to  Vickers  indented  glass  test  pieces  is 
plotted  in  Figure  6.  The  critical  applied  stress  am  defines  a  new  strength  level  [28] 

S  ■  Kc/2Yc^  ,  (6) 

which  is  always  less  than  the  evaluation  from  equation  (5).  Hence  the  presence  of  resi¬ 
dual  stresses  has  a  clearly  detrimental  effect  on  the  structural  integrity  of  materials. 

It  can  be  shown  from  the  fracture  mechanics  that  cm  «  (P/Kc)ys,  independent  of  the  initial 
flaw  size  Co,  so  that  control  of  strength  can  be  exercised  via  the  indentation  load  P. 
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It  is  this  feature  which  establishes  equation  (6)  as  an  attractive  base  formula  for  the 
accurate  determination  of  material  fracture  parameters  from  strength  data  [4] . 


Figure  6  -  Stable  Crook  Extension  in  Vickers-Indented  Glass  Prior 
to  Failure  in  Tensile  Loading.  Stress  and  oraok  size 
normalised  to  coordinate  values  at  maximum  in  curve. 
Closed  symbols  denote  specimens  confined  to  inert 
environment  prior  to  failure  test  to  suppress  poet- 
indentation  crack  creep ,  open  symbols  denote  corrosive 
(moist)  environment.  Note  failure  conditions  insensitive 
to  initial  crack  size.  Courtesy  of  D.B.  Marshall. 


One  area  of  particularly  intense  recent  interest  in  relation  to  the  issue  of 
materials  evaluation  is  that  of  "fatigue",  which  in  the  ceramics  testing  fraternity  is 
taken  to  mean  strength  loss  due  to  environmentally-assisted  crack  growth  during  a  static 
or  monotonically  increasing  applied  tensile  stress.  The  introduction  of  indentation 
flaws  into  fatigue  test  pieces  allows  for  determination  of  appropriate  rate  parameters 
for  the  underlying  crack  growth  process,  with  unprecedented  experimental  simplicity  and 
specimen  economy  [29,30].  More  significantly  in  the  present  context,  however,  it  high¬ 
lights  the  dramatic  increase  in  susceptibility  to  fatigue  effects  that  flaws  exhibit  when 
residual  stresses  are  present.  For  instance,  by  annealing  out  such  stresses  from  Vickers- 
indented  glass  test  pieces  the  lifetime  at  a  prescribed  applied  load  is  extended  by  some 
three  orders  of  magnitude  [31] .  There  are  important  implications  here  in  establishing 
design  criteria  for  brittle  materials,  especially  in  view  of  a  growing  conviction  that 
many  "natural"  flaws  nay  be  more  closely  simulated  by  indentation  cracks  in  their  "as- 
produced"  rather  than  "annealed"  state. 

At  the  time  of  writing,  studies  of  strength  degradation  for  indenters  subject 
to  translational  motion  have  only  just  begun.  Early  indications  are  that  the  linear 
flaws  produced  in  this  type  of  contact  are  likely  to  be  even  more  strongly  influenced 
by  residual  stress  effects  than  are  the  simple  point-contact  flaws  described  above. 


Erosive  Wear 


The  fracture  mechanics  relations  derived  earlier  can  also  be  used  to  develop  models  for 
erosive  wear,  where  one  seeks  to  inhibit  cracking,  and  for  machining,  where  one  seeks 
to  promote  cracking.  Sharp  indenters  are  most  effective  in  such  material  removal  pro¬ 
cesses  because  of  the  relative  ease  with  which  lateral  cracks  form  intersections  with 
the  free  surface.  Accordingly,  most  treatments  of  erosive  wear  and  machining  damage 
have  concentrated  on  this  contact  mode.  We  shall  consider  just  one  practical  example, 
that  of  erosion  by  normally  incident  sharp  particles,  to  illustrate  the  modelling 
procedure. 

The  basic  starting  equation  for  the  formulation  defines  the  approximate  amount 
of  material  removed  in  a  single  impact  event, 

AV  =  irc2d  ,  (7) 

where  c  is  the  radius  and  d  the  depth  of  the  subsurface  lateral  crack.  The  crack  radius 
is  estimated  from  the  lateral  crack  analogue  of  equation  (4) ,  and  the  depth  from  the 
penetration  of  the  deformation  zone.  An  appropriate  impulse  relation  is  once  again 
required  to  eliminate  contact  load  in  favour  of  incident  velocity.  Summation  over  all 
impacts,  assuming  zero  interaction  between  neighbours,  leads  to  a  functional  relation 
(usually  in  power- law  form)  for  the  total  volume  removal  rate  in  terms  of  the  particle 
kinetic  energy  and  material  parameters.  Unfortunately,  the  final  expression  obtained 
is  extremely  sensitive  to  minor  details  in  the  analysis  (note  equation  (7)  involves  the 
cube  of  a  linear  dimension,  whereas  the  strength  formulae  of  equations  (S)  and  (6)  in¬ 
volve  only  the  square  root),  and  erosive  data  are  notorious  for  their  scatter,  the 
result  of  which  is  an  accumulation  of  slightly  different,  experimentally  indistinguish¬ 
able  theories  [32,33].  Thus  the  complexity  of  the  erosion  process  admits  only  a  partial 
theoretical  account  of  experimental  behaviour. 
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Hardness  testing  is  gaining  popularity  as  a  scientific 
tool  for  evaluating  the  deformation  properties  of 
materials  (1],  Its  simplicity  and  economy  are 
particularly  well  suited  to  systematic  investigation 
of  the  many  variables  which  are  manifest  in  the 
general  deformation  response  of  any  given  material. 
Time  is  one  of  the  more  significant  of  these 
variables,  for  it  contains  the  key  to  the  kinetics 
and  dynamics  of  material  behaviour  at  the  micro¬ 
mechanical  level.  Such  rate  effects  could  be  an 
important  design  consideration  where  hardness 
values  are  taken  as  a  basis  for  predicting  strength, 
erosion  and  wear  characteristics  of  solids. 

Reports  in  the  literature  of  rate  dependencies 
in  hardness  measurements  are  sparse.  The  earliest 
studies,  on  metals,  showed  effects  which  appeared 
to  reflect  general  creep  properties;  thus  whereas 
indium  showed  substantial  decreases  in  hardness 
with  increasing  indentation  time  [2)  metals  with 
higher  melting  points  apparently  did  not,  unless 
the  test  temperature  was  raised  to  a  sufficiently 
high  level  [3].  Similar  indentation  creep  phenomena 
were  later  reported  for  some  ionic  solids  [4,  5|.  A 
study  on  a  wide  range  of  other  nonmetallic  solids 
showed  somewhat  analogous  behaviour  for  low- 
load  indentations  (<  IN)  [6],  but  the  effect 
disappeared  at  higher  loads;  this  led  the  authors 
of  that  study  to  conclude  that  they  were  observing 
a  spurious  manifestation  of  surface-adsorbed  water. 
Gunasekera  and  Holloway  [7],  in  perhaps  the  most 
detailed  investigation  of  this  kind,  found  strong 
rate  sensitivities  in  the  hardness  of  silicate  glass, 
depending  on  both  the  environment  and  the  state 
of  the  surface;  the  loads  used  in  their  study  were 
sufficiently  large  (1  to  5  N)  for  these  workers  to 
convince  themselves  they  were  not  simply  record¬ 
ing  an  artifact  of  either  the  experimental  procedure 
or  the  material. 

The  results  to  be  described  in  the  present  paper 
represent  preliminary  findings  of  a  proposed 


comprehensive  programme  to  investigate  time- 
dependent  mechanical  properties  of  selected 
engineering  materials.  Specifically,  we  report  on 
the  hardness  response  of  a  copper  single  crystal 
[(110)  surface  orientation],  polycrystalline  tung¬ 
sten  (grain  size  1  pm)  and  crown  silicate  glass 
(thermally  tempered,  to  minimize  the  incidence  of 
radial  cracking  at  the  indentation  corners);  these 
three  candidates  were  chosen  simply  to  embrace 
a  range  of  material  types. 

The  hardness  testing  apparatus  used  in  the 
study  represents  a  modification  of  that  described 
by  Gunasekera  and  Holloway  [7],  An  electro¬ 
magnetic  transducer  is  used  to  drive  the  indenter 
arm  vertically  downward  onto  the  surface  of  a 
rigidly  mounted  specimen.  The  compliance  of  the 
mechanical  support  for  the  coil  in  this  system 
is  sufficiently  small  that  the  load  delivered  is 
determined  exclusively  by  the  electrical  input  to 
the  transducer,  independent  of  any  indenter 
movement  relative  to  the  specimen  surface.  In 
our  experiments  a  sinusoidal  signal  is  used  to 
minimize  complications  due  to  inertial  effects.  The 
indenter  is  a  standard  Vickers  diamond  pyramid. 
The  indentation  load  is  monitored  directly  by 
means  of  a  piezoelectric  or  strain-gauge-instru¬ 
mented  cell,  depending  on  whether  the  contact 
period  is  less  than  or  greater  than  ^  1  sec.  These 
cells  are  generally  located  immediately  below  the 
specimen  (see  below,  however).  The  load-time 
pulses,  along  with  the  corresponding  activation 
pulses  from  the  sine-wave  generator,  are  recorded 
on  an  oscilloscope.  With  this  arrangement  an 
accuracy  of  2%  in  the  load  measurement  is 
readily  attainable.  The  indentations  produced 
are  then  measured  optically  in  the  usual  way,  and 
the  hardness  evaluated  from  the  averaged  impres¬ 
sion  diagonals  on  a  projected  contact  area  basis. 

Initial  runs  to  test  the  performance  of  the 
above  system  soon  revealed  difficulties  due  to 
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Figure  1  Vickers  hardness  as  function  of  contact  period 
for  three  materials.  Solid  lines  are  empirical  fits  to  the 
data. 

resonance  effects  at  contact  times  much  less  than 
10  msec.  In  this  region  the  load  pulse  began  to 
reduce  in  amplitude  and  to  lag  in  phase  relative  to 
the  input  signal;  spurious  oscillations  also  became 
apparent.  The  inertia  of  the  indenter  drive  system 
was  certainly  a  contributing  factor  here,  as  could 
be  demonstrated  by  removing  the  piezoelectric 
load  cell  from  below  the  specimen  and  relocating 
it  in  the  actual  indenter  arm;  the  recorded  peak 
load,  and  thence  the  hardness  evaluation,  were 
much  reduced  by  this  interchange.  We  may  note 
in  passing  that  such  artifically  low  readings  could, 
if  not  properly  accounted  for,  lead  one  to  conclude 
that  a  given  material  is  actually  undergoing  a 
rate-induced  softening,  e.g.,  due  to  localized 
adiabatic  heating  [8].  In  any  event,  we  avoid  such 
potential  ambiguities  in  interpretation  here  by 
including  only  data  from  indentations  whose  load 
pulses  were  entirely  free  of  any  resonance  effects. 

For  the  runs  proper,  indentations  were  made  on 
mirror-finished  surfaces  of  the  three  test  materials, 
in  air,  at  loads  which  could  be  considered  to 
sample  the  bulk  deformation  properties,  namely, 
*  5  N  for  the  metals  and  2.5  N  for  the  glass.  The 


results  are  shown  in  Fig.  1.  Each  data  point  in 
this  figure  represents  the  mean  and  standard 
deviation  of  measurements  from  five  impressions 
at  a  fixed  input  signal  setting.  Significant  increases 
in  hardness  are  observed  with  declining  contact 
period  for  the  three  months,  particularly  for  the 
tungsten  and  the  glass.  The  trend  for  the  latter  is 
commensurate  with  that  observed  by  Gunasekera 
and  Holloway  [7  ].  It  is  interesting  to  note  that  of 
the  two  metals  studied  it  is  the  one  with  the  lower 
melting  point,  copper,  which  is  relatively  insensitive 
to  time  variations. 

The  data  contained  in  Fig.  1  are  acknowledged 
to  be  restricted  in  several  respects:  the  number  of 
solids  investigated  is  hardly  large  enough  to  allow 
general  conclusions  to  be  drawn  concerning  the 
role  of  material  parameters  in  the  deformation 
kinetics;  the  range  of  contact  times  is  limited, 
such  that  any  extrapolations  into  the  realms  of 
either  dynamic  impact  or  long-term  loading  would 
carry  an  element  of  extreme  uncertainty;  and  no 
attempt  has  been  made  to  control  the  environ¬ 
ment.  Nevertheless,  the  results  serve  to  confirm 
that  hardness  can  indeed  be  subject  to  rate  effects. 
Furthermore,  systematic  study  of  such  trends 
might  be  expected  to  throw  some  light  on  the 
micromechanics  of  the  processes  which  control 
the  deformation  properties  of  materials. 
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A  theory  is  formulated  for  the  general  fatigue  response  of 
brittle  flaws  which  experience  residual  stress  concentrations. 
The  indentation  crack  is  taken  as  a  model  flaw  system  for  the 
purpose  of  setting  up  the  basic  fracture  mechanics  equations, 
but  the  essential  results  are  expected  to  have  a  wider  range 
of  applicability  in  the  strength  characterization  of  ceramics. 
A  starting  fatigue  differential  equation  is  first  set  up  by 
combining  an  appropriate  stress  intensity  factor  for  point-  or 
line-contact  flaws  with  a  power-law  crack  velocity  function. 
Analytical  solutions  are  then  obtained  for  the  case  of  static 
fatigue.  The  resulting  relation  between  lifetime  and  failure 
stress  is  shown  to  have  exactly  the  same  power-law  form  as  the 
conventional  solution  for  Griffith  (residual-stress-free)  flaws. 
This  “equivalence”  is  used  as  a  basis  for  extending  the  results 
to  dynamic  fatigue.  A  comparison  of  these  analytical  solutions 
with  numerical  counterparts  defines  the  limits  of  accuracy  of 
the  theoretical  procedure.  However,  while  the  form  of  the  life¬ 
time  relation  remains  invariant,  the  values  of  the  exponent  and 
coefficient  differ  significantly  for  flaws  with  and  without  re¬ 
sidual  stress.  Accordingly,  the  application  of  conventional  fa¬ 
tigue  theory  to  evaluate  crack  velocity  parameters,  without  due 
regard  for  the  nature  of  the  critical  flaw,  can  lead  to  serious 
errors.  Explicit  conversion  formulas  are  given  for  trans¬ 
forming  “apparent”  velocity  parameters  for  indentation  flaws 
directly  into  “true”  parameters.  The  implications  of  these 
results  concerning  the  use  of  the  indentation  method  for  mate¬ 
rials  evaluation  are  discussed. 

I.  Introduction 

IT  is  well  recognized  that  the  failure  of  brittle  materials  is  gov¬ 
erned  by  the  micromechanics  of  crack  growth  from  small  flaws, 
and  that  chemical  enhancement  of  this  crack  growth  can  cause 
significant  reductions  in  the  strength  with  increasing  time  under 
load.  Embodied  in  the  conventional  fracture  mechanics  approach 
to  “fatigue"  phenomena  of  this  kind' 3  are  three  underlying  as¬ 
sumptions:  (a)  The  time  dependence  of  the  loading  stresses,  taken 
to  act  uniformly  across  the  prospective  crack  plane,  is  specifiable; 
( b )  the  driving  force  on  the  extending  crack  is  uniquely  determined 
at  any  given  characteristic  length  by  these  applied  loading  stresses; 
(c)  the  rate  of  crack  extension  is  in  tum  uniquely  determined  by 
some  well-defined  function  of  the  driving  force  for  any  given 
material/environmcnt  system.  These  assumptions  allow  one  to 
write  down  a  differential  equation  in  crack  length  and  time,  the 
solution  of  which  defines  the  stress  conditions  at  failure.  The 
widespread  success  enjoyed  by  the  fracture  mechanics  formulation 
arises  from  the  amenability  to  solutions  in  simple,  closed  form, 
which  provide  a  convenient  basis  for  lifetime  predictions. 

Apart  from  the  clear-cut  distinction  made  between  loading  at 
constant  stress  (“static  fatigue”)  and  constant  stress  rate  ("dynamic 
fatigue”),  surprisingly  little  attention  has  been  devoted  to  the  ef¬ 
fects  that  potential  variations  in  the  starting  equations  may  have  in 
the  lifetime  analysis.  Wiederhom  and  Ritter1  examined  the  crack 
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velocity  function,  and  concluded  that  any  of  the  commonly  ac¬ 
cepted  empirical  forms  may  fit  fatigue  data  equally  well  (although 
extrapolations  beyond  the  data  range  could  lead  to  significant  dis¬ 
crepancies  in  the  predictions).  Of  the  assumptions  listed  above,  the 
second  has  been  subjected  to  least  scrutiny,  it  generally  being 
assumed,  without  question,  that  the  flaws  respond  in  the  classical 
“Griffith"  sense;  that  is,  the  flaw  is  driven  solely  by  the  applied 
loading,  this  force  increasing  monotonically  with  the  crack  size 
until  some  instability  condition  is  met.4'1  It  is  thus  implicit  in  the 
statement  of  the  problem  that  any  preexisting  stresses  which  may 
have  been  responsible  for  generating  the  critical  flaw  in  the  first 
place1'6  have  long  since  ceased  to  be  a  significant  contributing 
factor  in  the  net  driving  force  on  the  system. 

However,  recent  studies  of  controlled  flaws  produced  by  inden¬ 
tation  in  strength  test  pieces  have  demonstrated  that  residual  crack - 
generation  stresses  can  have  a  profound  influence  on  the  crack 
evolution  to  failure.7  The  source  of  the  residual  field  in  this  case 
is  elastic-plastic  mismatch  at  the  boundary  of  the  deformation  zone 
which  encases  the  sharp  point  and  edges  of  the  contacting  body.” 
It  then  becomes  necessary  to  incorporate  a  residual-contact  term 
into  the  fracture  mechanics  equation  for  the  crack  driving  force. 
Characteristically,  this  contribution  decreases  monotonically  with 
crack  size.”  The  resulting  expression  for  the  net  force  on  the  crack 
now  takes  on  a  considerably  more  complicated  form.  This  com¬ 
plexity  is  such  that  the  appropriate  fatigue  differential  equation  no 
longer  appears  to  have  simple  analytic  solutions.  Accordingly,  the 
first  systematic  investigations  of  residual-stress  effects  in  fatigue, 
using  results  from  dynamic9  and  static10  loading  tests  on  Vickers- 
indented  soda-lime  glass  in  water  as  a  data  base,  were  made  by 
obtaining  numerical  solutions  specific  to  one  particular 
indentcr/material/environment  system.  A  subsequent  analysis," 
based  on  a  reformulation  of  the  differential  equation  in  terms  of 
judiciously  normalized  variables,  allowed  for  generalization  of  the 
numerical  procedure  to  include  solutions  for  all  possible  systems. 
Most  notably,  this  last  study  produced  an  empirical  dynamic  fa¬ 
tigue  relation,  for  flaws  satisfying  a  power-law  crack  velocity 
function,  which  was  indistinguishable  in  form  from  that  derived 
analytically  for  Griffith  flaws.  The  exponents  and  coefficients  in 
this  relation  were  not,  however,  identical  in  the  two  cases;  in 
particular,  the  values  of  the  fatigue  exponent,  which  for  Griffith 
flaws  is  a  direct  measure  of  the  corresponding  exponent  in  the 
crack  velocity  function,  differed  by  some  30%.  A  case  study  on  a 
glass-ceramic13  confirmed  these  and  other  features  of  the  residual- 
stress  theory,  and  outlined  several  unique  advantages  of  the  inden¬ 
tation  method  as  a  means  for  evaluating  basic  fatigue  parameters. 

One  point  that  must  be  made  at  the  outset  is  that  indentation 
cracks  should  not  be  regarded  simply  as  artificially  introduced 
entities  which  bear  no  resemblance  to  strength-controlling  flaws  in 
real  materials.  There  is  growing  evidence  that  the  degrading  sur¬ 
face  damage  which  many  ceramic  components  experience  in  fin¬ 
ishing  (e.g.  machining)11  or  in  service  (sharp  paiticle  impact)14  are 
characterized  by  the  same  residual  stress  effects  as  are  indentation 
flaws.  Indeed,  the  observation  of  strongly  analogous  local  stress 
field  effects  about  microstructural  flaws  in  ceramics”  suggests  that 
the  presence  of  residual  crack  driving  forces  may  be  the  rule  rather 
than  the  exception. 

Viewed  against  this  background  the  solution  of  the  indentation 
fatigue  problem  takes  on  a  broader  significance.  Accordingly,  the 
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Fi«.  I.  Indentation  flaw;  crack  of  size  c  is  formed  at 
contact  load  P,  and  subsequently  subjected  to  applied 
tensile  stress  < r„.  Deformation  zone  about  contact  gives 
rise  to  residual  stress  field  which  contributes  tn  crack 
driving  force. 


numerical  base  of  the  earlier  analyses'*- "  must  be  seen  as  re¬ 
strictive.  ideally,  one  would  like  to  be  able  to  obtain  analytical 
solutions  of  the  master  differential  equation  in  its  most  general 
possible  form.  Such  solutions  would  provide  a  sounder  basis  for 
making  intercomparisons  between  (a)  flaws  of  different  geometry, 
e.g.  “point”  flaws  produced  in  normal  loading  vs  “line”*  flaws 
produced  in  sliding  loading;  (b)  different  crack  velocity  functions; 
(c)  static  vs  dynamic  fatigue.  In  this  paper  we  present  an  analysis 
which  meets  this  ideal  at  least  in  part,  the  greatest  restrictions  being 
the  need  to  retain  a  power-law  crack  velocity  function  and  to  obtain 
the  dynamic  fatigue  solutions  by  an  “equivalence"  argument.  The 
ensuing  fatigue  relation  between  lifetime  and  failure  stress  con¬ 
firms  the  findings  of  the  previous  empirical  studies,  but  now 
provides  more  explicit  expressions  for  obtaining  crack  velocity 
parameters  from  fatigue  plot  slopes  and  intercepts. 

II.  Analytical  Solution  of  Fatigue  Differential  Equation  for 
Constant  Applied  Stress 

(!)  The  Stress  Intensity  Factor  and  the  Inert  Strength 
The  key  step  in  generalizing  the  Griffith-flaw  concept  to  include 
residual-stress  effects  is  an  appropriate  expression  for  the  crack 
driving  force.  The  essential  variables  expected  to  appear  in  any 
such  expression  are  depicted  in  Fig.  I :  c  is  the  characteristic  crack 
size,  P r  the  indentation  load  which  determines  the  level  of  the 
residual  field,  and  o„  the  subsequently  applied  tensile  stress  which 
takes  the  system  to  failure.  In  this  work  we  shall  make  a  special 
distinction  between  the  geometrical  extremes  of  “point”  and  “line” 
flaws;  the  former  defines  a  crack  configuration  of  semicircular 
profile  centered  about  a  point-force  contact  (F,  “force),  the  latter 
a  crack  of  straight  front  parallel  to  a  line-force  contact  (F,“force/ 
length).  We  shall  also  assume  that  the  cracks  are  "well  developed,” 
i.e.  are  large  compared  to  the  deformation  zone  from  which  these 
cracks  initiate.  Following  our  previous  procedure,"  the  requisite 
crack  driving  force  may  be  formulated  in  terms  of  the  composite 
stress  intensity  factor 

K=Kr+K.  (la) 

where  the  terms 


Fig.  2.  Stress  intensity  factor  as  function  of  crack  size 
(Eq.  (I))  for  point  (r= 3)  and  line  (r=l)  flaws.  Inclined 
broken  lines  represent  individual  residual  and  applied 
components  (Eqs.  (lb)  and  (lr)):  solid  curves  represent 
composite  functions. 


minimum,  tending  asymptotically  on  either  side  of  this  minimum 
to  plots  of  the  residual  component  of  the  stress  intensity  factor  in 
the  small-crack  limit  and  of  the  applied  component  in  the  large- 
crack  limit.  The  coordinate  variables  are  normalized  to  the  refer¬ 
ence  points  M,  which  are  of  special  significance  in  establishing 
baseline  levels  for  the  fatigue  strength  characterization";  in  addi¬ 
tion,  this  mode  of  plotting  foreshadows  the  normalization  scheme 
to  be  adopted  in  a  later  section  of  the  paper.  For  arbitrary  values 
of  applied  stress  the  condition  dK/dc= 0  defines  the  minimum  in 
Eq.  (I),  which  we  designate  by  asterisk  notation: 


K* = (r+1  )x,F,/c*2 = [( r + 1  )/r]dwr« c'f 

(2a) 

C*  =  (rXrPr/<Pr<T.f^" 

(2b) 

The  curves  with  their  minimum  at  M  correspond  to  the  special  case 
K*=Kr\  at  this  point  we  may  appropriately  identify  the  critical 
variables  oa-crm,  c*=c„,  i.e. 

<r„=[r/(r+ 1  )'** ""]*T "7  tM*/*)"' 

(3  a) 

cm=[(r+\)x,P./K,T' 

(3 b) 

It  may  be  shown  from  Eqs.  (2)  and  (3)  that 

Kn/K'-lvJotjr*" 

(4 a) 

c*/  cm—(<rm/ trjf'r* 1  ’ 

(4b) 

K'-X'P'/c”2  (lb) 

K.=4 Ikhcm  GO 

respectively  represent  the  contributions  from  the  residual  contact 
field  and  the  applied  loading:  here  r=3  for  point  flaws  and  r=  1  for 
line  flaws;  and  4>,  are  dimensionless  parameters  of  the  inden¬ 
tation  stress  field  and  the  crack  geometry,  respectively.  This  for¬ 
mulation  is  subject  to  some  modification  due  to  the  influence  of 
such  factors  as  secondary  crack  systems,  spurious  surface  stress 
states,  etc.;  detailed  discussion  of  these  factors  is  available  else¬ 
where,"  and  will  consequently  be  pursued  no  further  here. 

Equation  (I)  has  certain  features  which  are  pertinent  to  the 
fatigue  analysis  to  follow.  These  features  arc  evident  in  the  plots 
of  the  function  K(c)  in  Fig.  2  for  both  point  and  line  flaws.  The 
curves  for  fixed  values  of  P,  and  a,  pass  through  a  well-defined 


so  that  as  a r„  drops  below  om  the  position  of  the  minimum  in  Fig.  2 
displaces  downward  and  to  the  right  relative  to  M 
The  points  at  which  the  curves  in  Fig.  2  intersect  the  horizontal 
line  K=KC  correspond  to  equilibrium  crack  configurations,  stable 
or  unstable  according  to  whether  the  branches  have  negative  or 
positive  slope.  The  stable  equilibria  define  appropriate  initial  con¬ 
ditions  for  ensuing  fatigue  fracture;  c0  for  loading  at  constant  stress 
rate,  Ci  at  constant  stress.  In  practice,  fatigue  effects  will  be  mani¬ 
fest  in  the  postindentation  crack  configuration  before  application  of 
tensile  loading,  causing  subcritical  extension  from  cn  to  some  non¬ 
equilibrium  size  Co  (which  may  or  may  not  exceed  o).  The  unstable 
equilibria  at  c,  define  the  final  configuration  immediately  prior  to 
the  onset  of  catastrophic  failure.  It  is  useful  at  this  stage  to  point 
out  that  an  instability  configuration  can  be  achieved  without  ever 
departing  from  an  equilibrium  state,  by  steadily  increasing  the 
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NotmaUzad  Failure  Sins*,  S, 

Fig.  3.  Normalized  plots  of  lifetime  vs  failure  stress  for  point 
and  line  flaws;  static  fatigue  results.  Data  points  and  solid  lines 
represent  numerical  and  analytical  solutions  of  master  fatigue 
differential  equation,  respectively.  Note  inert  strength  cutoff  at 


Nonwfcfd  Fathire  Streaa.  8, 

.  4.  Normalized  plots  of  lifetime  vs  failure  stress  for  point 
line  flaws;  dynamic  fatigue  results.  Data  points  and  solid 
lines  represent  numerical  and  analytical  solutions,  respectively. 


applied  stress  <x„  from  zero  to  <r„,  thereby  causing  the  crack  to 
grow  stably  from  c0  to  c„;  this  takes  us  along  K=KC  to  point  M  in 
Fig.  2,  where  c,  merges  with  c,  to  produce  spontaneous  failure. 
Since  such  conditions  are  most  closely  met  in  nonreactive  test 
environments  we  may  define  an  “inert”  strength  ov=tr,,  which 
from  Eq.  (3)  can  be  written  in  the  form 

<T.=[r/(r+\)]K./<l>,c”  (5) 

independent  of  initial  crack  size. 

For  comparison,  the  corresponding  inert  strength  a?  for  Griffith 
flaws  follows  from  Eq.  ( I )  in  the  limit  of  *,=0,  with  spontaneous 
failure  occurring,  without  any  precursor  crack  growth,  at  ch  say: 

a?=Kj4>rc!'2  (6) 

In  this  case  the  strength  is,  of  course,  sensitive  to  the  initial  flaw 
size;  the  appropriate  value  of  c„  i.e.  c0,  ch,  or  c(  (or  indeed  any 
other  such  crack  dimension),  will  depend  on  the  nature  of  the 
mechanical,  thermal,  and/or  chemical  (or  other)  processes  re¬ 
sponsible  for  removing  the  residual  contact  field  between  crack 
formation  and  strength  testing. 

(2)  Formulation  and  Solution  of  the  Fatigue  Equation 
The  fracture  mechanics  approach  to  the  fatigue  problem  begins 
with  the  assumption  that  a  crack  velocity  function  may  be  written, 
for  a  given  material/environment  system,  in  the  form  v=v{K).  In 
combination  with  die  stress  intensity  factor,  K=K(r,  P„  <r«,c),  in 
Eq.  ( I )  and  the  specified  time  variation  of  the  applied  tensile  field, 
ir„=cr„(f),  the  velocity  function  assumes  the  form  of  a  differential 
equation,  dc/dt- v(r,  P„  <r.(t),  c].  This  equation  must  be  solved  for 
the  time-to-failure,  //,  needed  to  take  the  crack  from  c,  to  cf,  at 
which  point  the  stress  level  defines  the  fatigue  strength,  er.^oy. 
The  primary  objective  of  any  such  analysis  is  to  determine  the 
relation  between  trf  and  tf  (or  some  equivalent  parameter,  such  as 
stress  rate  er,  in  dynamic  fatigue). 

To  proceed  with  the  solution  of  the  differential  equation  it  is 
necessary,  of  course,  to  know  the  form  of  the  crack  velocity  func¬ 
tion.  In  this  paper  wc  adopt  the  simple  power-law  relation 

v-vo (K/K'T 


where  v0  and  n  are  empirical  quantities.  At  the  moment,  this  is  the 
only  one  of  the  commonly  used  crack  velocity  functions  for  which 
we  have  been  able  to  obtain  analytical  solutions  for  flaws  with 
residual  stress.  Thus  we  obtain 

dc/dt=(va/KM[r.  Pr.  cr.(f).  c]}'  (8) 

as  our  master  starting  equation. 

(A)  Griffith  Flaws;  Static  and  Dynamic  Fatigue  Solutions: 
The  standard  solutions  of  Eq.  (8)  for  Griffith  flaws  are 
well  known.1"3  It  is  nevertheless  instructive  to  include  them  here 
as  a  base  for  comparing  later  solutions.  Inserting  *,=0  into 
Eq.  (1)  yields 

dc/dt=(v0/K7)K2 

=(v0 fKMurAW  (9) 

The  problem  accordingly  reduces  to  one  of  straightforward  inte¬ 

gration  by  separation  of  variables.  For  static  fatigue,  i.e. 
cr.= cons  tan  1=0/,  Eq.  (9)  becomes  simply 

Jdt*[tf/v dele*2  (10) 

The  integrated  result  may  be  written  in  the  form 

t/of-A,  (II) 

where,  in  the  usual  approximation  (c#/<y)<""2*'2  <  I ,  and  in  conjunc¬ 
tion  with  Eq.  (6),  we  obtain 

A,  *  2K’/(n  -  2)v0d'?c/"“ 21,2 

=[2/(n-2)]oi°" c(/vo  (12) 

The  appearance  of  c,  as  the  controlling  crack  size  in  Eq.  (12) 

reflects  the  fact  that  the  crack  system  spends  the  greater  propor¬ 
tion  of  its  time  evolution  to  failure  in  the  region  of  minimum 
driving  force. 

The  corresponding  solution  for  dynamic  fatigue,1"'  i.e.  for 
with  constant  can  be  expressed  by  the  same  power  law 
as  Eq.  (II),  but  with  A,  replaced  by  A*  where 

(13) 


<7>  ax  A,-(n+l)A. 
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IB)  Flaws  with  Residual  Stress;  Static  Fatigue  Solution:  For 
Che  general  case  where  fO  in  Eq.  (1),  Eq  (8)  becomes 

de/dt — ( v«/ AT )[  A'.  <•/£■„]” 

=(vlt/K:ix,Pr/c"2  +  'i>,<TM)c'2]''  ( 14) 

Iniegration  by  separation  of  variables  is  no  longer  straightforward 
as  it  was  for  Griffith  Haws.  Even  in  static  fatigue,  for  which  the 
right  side  of  Eq.  (14)  contains  no  explicit  terms  in  time,  it  is  not 
immediately  clear  how  the  failure  stress  might  be  extracted  from 
the  integral  and  thereby  related  to  lifetime.  It  was  difficulties  of  this 
kind  which  motivated  the  numerical  approach  described  in 
Ref.  II. 

However,  it  can  now  be  demonstrated  that  Eq.  ( 14)  does  have  an 
analytical  solution  at  (r„=constant  =  o/.  Since  the  integrated  life¬ 
time  is  controlled  at  one  extreme  by  K,  for  small  cracks  (c<c»)  and 
at  the  other  extreme  by  K„  for  large  cracks  (cS>c*).  it  is  appropriate 
to  introduce  a  variable  that  defines  the  relative  stress  intensity 
factor  at  either  extreme.  We  accordingly  choose 

(=KJ(K,+KJ 

=  l/[l  HXrPrhh.cr,)/^"1]  (15) 

where  £  defines  the  fraction  of  the  total  stress  intensity  factor  wriich 
is  associated  with  the  applied  field,  this  fraction  increasing  mono- 
tonicaJly  with  crack  size.  A  second,  important  motivation  for 
introducing  this  reduced  variable  is  to  convert  the  lifetime  integral 
to  a  dimensionless  form  so  as  to  display  the  applied  stress  and 
indentation  load  dependence  separately.  After  considerable 
manipulation.  Eq.  (14)  becomes 


r- 

j  0 


{I2*r/(r+ 1  )*]/[(*<*)•' 


X  £"•**»’>»  ,(l-f)'"‘2"~,,“,df  (16) 

and  we  note  that  try  does  indeed  appear  outside  the  integral .  At  this 
stage  a  significant  simplification  in  the  expression  can  be  made  by 
introducing  the  quantity 

n’=(rn+2)/iri- 1 )  (17) 

This  will  be  seen  later,  with  the  benefit  of  hindsight,  to  be  a 
particularly  convenient  choice  of  substitution.  A  further  sim¬ 
plification  is  to  assign  the  following  values  to  the  limits  of  integra¬ 
tion  in  Eq.  (16),  consistent  with  our  identification  of  the  initial  and 
final  conditions  in  Fig.  2  and  Eq.  (IS), 

(i=il/rO/cl*/K,—*0  (18a) 

(f=  (186) 

in  which  case  the  integral  in  Eq.  (16)  reduces  to  the  beta  function.* 
These  value  assignments  are  tantamount  to  saying  that  the  initial 
conditions  are  governed  by  the  residual  component  of  the  stress 
intensity  factor,  the  final  conditions  likewise  by  the  applied  compo¬ 
nent.  In  any  case,  the  final  result  is  not  expected  to  be  sensitive  to 
the  limits  of  integration,  bearing  in  mind  that  the  region  of  mini¬ 
mum  driving  fone  in  Fig.  2,  which  must  control  the  fracture 
kinetics,  lies  intermediate  to  the  initial  and  final  configurations. 
This  assumption  is  verified  in  the  Appendix.  With  these  sim¬ 
plifications  Eq.  (16)  reduces  to 

<r=[2K,7 (r  *  I  »„<«k,<r,)" (*,P,)"  *'  ]B(n\ n-n')  ( 19) 

with  the  beta  function 


B(n\ «-/»•)=  [  f -'(!-£)" 

•'(i 


Thus  the  requisite  tciation  between  lifetime  and  failure  stress  has 
the  familiar  form 

1/0,"=*:  (2D 


•See  any  xandard  km  no  advanced  mathematical  methods 


where  we  obtain,  in  conjunction  with  Eqv  (t)  and  (5), 

A.'  •  i  2F"  »  1 },-.  .!//:<*  P.\"  "  ]/((«',  n-  n') 

~  i!  2tr  ■  If  •;•>"  \lhn\  n-n'  )}it"  i  „,/r„ 

~[Kjr/tr-  I  )/i'  j ' (22) 

the  Iasi  expression  arising  from  an  asymptotic  expansion  of  the  beta 
function  lor  laree  n  (Appendix). 

While  the  form  of  bq.  (21)  is  indistinguishable  from  that  of 
Eq.  (II)  for  Griffith  flaws,  the  exponent  and  coefficient  are  sig¬ 
nificantly  modified.  The  implications  of  these  modifications  are 
discussed  in  Section  IV 

III.  Comparison  of  Analytical  and  Numerical  Solutions: 

Extension  to  Constant  Stress-Rate  Loading 

It  is  instructive  to  examine  the  mutual  consistency  of  the  results 
obtained  analytically  in  Section  II  and  numerically  in  Ref.  II.  To 
do  this,  and  to  broaden  the  scope  of  the  treatment  to  date,  numeri¬ 
cal  solutions,  previously  confined  to  the  special  case  cr„=constant 
and  r=3."  are  here  generated  for  each  combination  of  fatigue 
(static  or  dynamic)  and  flaw  (point  or  line)  types. 

Following  Ref.  1 1 ,  reduced  variables  are  introduced  as  follows: 


Su-<rja„  (23a) 

C=c/cm  (236) 

T=tvjcm  (23c) 

Reference  to  Eq.  (5)  then  allows  the  differential  equation  for  flaws 

with  residual  stress.  Eq.  (14),  to  be  expressed  in  a  more  universal 
form, 

dC/dT=  { 1  /(r+  1  )C"2  +[r/(r+ 1  )]S.(DC''2}"  (24) 


which  is  especially  amenable  to  numerical  analysis.  This  equation 
is  solved  by  a  stepwise  integration  procedure  for  the  (reduced) 
time-to-failure  7)  to  take  the  crack  from  C,  (approximated  by  the 
value  C0=\/(r+\fr  at  S„=0)  to  C>."  The  corresponding  failure 
stress  Sf  at  the  critical  end  point  of  the  integration  is  seen  to  be 
uniquely  determined  by  the  values  of  r  and  n  and  the  form  of  the 
stressing  function  SJT).  It  is  noted  that  the  normalization  scheme 
conveniently  relates  all  variables  to  the  inert  strength  state,  i.e. 
S,-S„=  1  and  C„=  1 ,  as  represented  by  the  reference  point  M  in 
Fig.  2. 

To  facilitate  the  required  comparisons,  let  us  translate  the  essen¬ 
tial  results  of  the  analytical  treatment  in  Section  II  into  reduced 
variable  notation.  Taking  static  fatigue  first,  i.e.  S„=constant=S/t 
Eq.  (21)  retains  its  basic  form. 

W=A;  (25) 

with  the  coefficient  relating  to  its  counterpart  in  Eq.  (22)  as 

A,'  =  A.'vn/ tX?Cm 

=[2(r+l)""'/r"]B(n,  n-n') 

a[8rr/<r+l)]’7n"'2  (26) 

where  the  last  line  in  this  equation  is  an  asymptotic  approximation 
in  n'  (Appendix). 

At  this  point,  it  is  pertinent  to  recall  the  essential  equivalence  of 
the  static  fatigue  solutions  for  flaws  with  and  without  residual 
stress,  Eqs.  (21 )  and  ( 1 1 ).  It  would  appear  reasonable  to  expect  this 
equivalence  to  extend  to  dynamic  fatigue,  Sa-S„T  (S„  constant). 
Accordingly,  since  in  the  case  of  Griffith  flaws  the  static  and 
dynamic  solutions  can  be  shown  analytically  to  have  exponents 
which  are  identical  and  coefficients  which  relate  linearly  via 
Eq.  (13),  we  may  proceed  by  analogy  and  modify  Eq.  (25),  thus 

W-AJ  (27) 

where  we  have  from  hq.  (26) 

A; -•  Aji-n/cr; c„-(n'  +■ 1  )A.' *[8ir/(r  +  I ))' V '  2  (28) 

with  the  same  asymptotic  approximation  as  made  in  Eq.  (26). 

Figures  3  and  4  show  the  results  of  both  analytical  and  numerical 
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Fig.  5.  Plots  showing  variation  of  exponent  n'  («,'  and  nj  cannot  be 
distinguished  in  these  plots)  and  coefficients  A,‘  and  A, I  with  n.  for  point 
and  line  flaws.  Data  points  and  solid  lines  represent  numerical  and  analyti¬ 
cal  evaluations,  respectively. 


calculations  on  logarithmic  plots  of  lifetime  vs  failure  stress  for 
static  and  dynamic  fatigue,  respectively,  and  for  point  and  line 
flaws.  The  data  points  in  these  plots  are  numerical  evaluations  of 
Eq.  (24)  at  either  fixed  S„  (Fig  3)  ori„  (Fig.  4)  for  selected  values 
of  n.  The  corresponding  solid  lines  are  analytical  representations 
of  Eqs.  < 25)  and  (27).  in  conjunction  with  Eq.  (17).  The  degree  of 
correlation  between  data  points  and  solid  curves  reflects  the  accu¬ 
racy  of  the  theory  developed  in  Section  II.  In  this  context  the 
tendency  for  systematic  departures  to  increase  at  low  n  (or  n') 
values  may  be  taken  as  a  measure  of  the  range  of  validity  of  the 
asymptotic  beta  function  expansion  used  in  the  evaluation  of  the 
integral  in  Eq.  ( 16).  Special  note  may  be  made  of  the  fact  that  the 
correlations  appear  to  be  as  strong  for  the  dynamic  fatigue  results 
in  Fig.  4  as  for  the  static  fatigue  in  Fig.  3,  thereby  providing  some 
justification  for  the  equivalence  argument  adopted  to  extend  the 
theoretical  analysis  cailier  in  this  section. 

A  more  detailed  comparison  of  results  is  obtained  by  considering 
the  specific  variations  of  the  exponent  n'  and  coefficients  A’  in 
Eqs.  (25)  and  (27)  with  n.  This  is  done  in  Fig.  5  for  point  and  line 
flaws.  The  data  points  again  represent  numerical  evaluations,  of 
Eq.  (24).  obtained  from  slope  and  intercept  evaluations  in  the 
"linear"  region  ol  the  lifetime  vs  failure  stress  plots;  the  solid 
curves  likewise  represent  analytical  evaluations  from  Eqs.  (17), 

(26),  and  (28).  The  disagreement  between  the  two  approaches 
is  <  1 '/<•  for  the  export  :nt  and  <I09£  for  the  coefficients  in  the 
domain  n>l(). 


IV.  Discussion 

The  formulation  presented  in  this  paper  demonstrates  that  analy¬ 
tical  solutions  can  be  ootained  to  the  fracture  mechanics  fatigue 
equations  for  flaws  which  are  subject  to  residual -contact  driving 
forces  and  which  extend  according  to  a  power-law  crack  velocity 
function  A  major  feature  of  these  solutions  is  the  fact  that  they  are 
identical  in  basic  form  to  those  obtained  for  Griffith  flaws  in 
conventional  fatigue  analysis.  Thus,  from  a  standard  linear  plot  of 
lifetime  vs  failure  stress  in  logarithmic  coordinates  it  would  not  be 
possible  to  determine,  without  independent  information  on  the 
crack  velocity  parameters,  whether  the  flaws  in  a  given  materia!/ 
environment  system  are  or  are  not  influenced  by  residual  stresses. 
This  conclusion  should  provide  some  comfort  to  those  who  have 
advocated  the  exclusive  use  of  strength  data  for  predicting  compo¬ 
nent  lifetimes,  since  the  nature  of  the  flaw  now  docs  not  enter  into 
consideration  (unless  perhaps  the  predictions  require  extrapolation 
beyond  the  data  range"’).  Of  particular  importance  in  this  context 
is  the  widely  proposed  use  of  dynamic  fatigue  testing  as  the  source 
of  base  data  for  design  requirements;  in  terms  of  a  lifetime  vs 
failure  stress  diagram  the  curve  for  static  loading  generates  from 
that  for  dynamic  loading  via  a  simple  connecting  relation  for  the 


intercepts.  To  illustrate.  Fig.  6  show  s  data  from  both  dynamic"  and 
static"'  fatigue  test  runs  on  indented  soda-lune  glass  in  water  A 
least-squares  fit  through  the  dynamic  fatigue  data  gives  n' .  Eq 

(27) ;  subsequent  translation  downward  through  logtii'  +  l).  Eq. 

(28) .  generates  the  static  fatigue  curve.  The  agreement  between 
prediction  and  observation  in  the  latter  case,  notwithstanding  the 
experimental  scatter,  may  be  taken  as  a  justification  for  the 
"equivalence"  argument  used  to  infer  the  existence  of  Eqs.  (27)  and 
(28)  from  their  Griffith  counterparts.  Eqs.  (II)  and  (13). 

There  is.  however,  considerable  danger  in  analyzing  fatigue 
results  from  materials  whose  flaw  characteristics  are  not  known. 
This  is  particularly  so  if  one  attempts  to  relate  the  strength  data  to 
the  parameters  in  the  crack  velocity  function.  The  point  is  most 
readily  demonstrated  by  considering  the  exponent  in  this  function 
in  terms  of  the  slope  of  lifetime  vs  failure  stress  plot.  From  Eq.  ( 1 1 ) 
for  Griffith  flaws  the  slope  gives  n  (more  strictly,  its  negative) 
directly,  whereas  from  Eq.  (21)  for  flaws  with  residual  stress  the 
corresponding  slope  is  n'\  in  the  latter  case  evaluation  of  the  true 
exponent  requires  inversion  of  Eq.  (17): 

«=4n'/3-2/3  (r=  3)  (29  a) 

n=2n’-2  (r=l)  (29b) 


Thus  substantial  errors  may  be  incurred  if  the  conventional  theory 
of  fatigue  is  used  to  analyze  data  for  contact-induced  flaws  which 
have  undergone  no  subsequent  stress  relaxation.  For  example,  in 
Fig.  6  the  apparent  exponent  is  determined  at  n'  =  13. 7^.0. 2, 
whence  Eq.  (29a)  appropriate  to  point  (Vickers-induced)  flaws 
predicts  n-  17  6±0.3;  this  latter  value  is  close  to  the  true  exponent 
n=I7.9±0.5  obtained  from  comparative  tests  on  indented  speci¬ 
mens  subjected  to  an  anneal  treatment  prior  to  strength  testing,* 
and  lies  in  the  range  of  I6<n  <  19  generally  found  for  large-scale 
cracks  in  the  system  soda-lime-glass/water. 16  The  discrepancy  be¬ 
tween  apparent  and  true  exponents  is  predicted  to  be  even  greater 
for  linear  flaws,  almost  a  factor  of  2,  Eq.  (29 b).  It  is  interesting  to 
note  that  discrepancies  of  this  order  were  found  by  Pletka 
and  Wiederhom1,  in  certain  ceramics  where  machining  damage 
provided  the  strength-controlling  flaws. 

Similar  care  must  be  exercised  in  evaluating  the  coefficient  of 
the  crack  velocity  function  from  the  intercept  on  a  fatigue  plot.  For 
flaws  with  residual  stress,  inversion  of  Eqs.  (26)  and  (28)  gives 


v(,=  A  'rr"cm/\'  (30) 

where 

A'  =[2ir/n')<17  (tr„ = constant,  r= 3)  (31a) 

A’=(47r/n’)1,J  ( <ja- constant.  r=  1)  (31b) 

=(2im')'11  (cr„=constant,r=3)  (31c) 

A'=(4irn')''2  (cr„=constant,r=l)  (3ld) 


Inversion  of  the  corresponding  intercept  relations  for  Griffith 
flaws,  Eqs.  (12)  and  (13),  leads  to  an  analogous  result,  but  with 
a  different  A'  term  and  with  c,  replacing  c„  as  the  controlling 
crack  parameter. 

We  should  emphasize  here  that  Eq.  (31)  represents  only  a  first- 
order  approximation;  additional  terms  in  a  series  expansion  may  be 
required  where  accuracy  requirements  are  stringent  (Appendix). 
However,  since  lifetimes  are  generally  plotted  in  logarithmic  coor¬ 
dinates  the  present  approximation  will  usually  suffice,  except  per¬ 
haps  at  nSlO. 

This  sensitivity  of  the  slope  and  intercept  terms  in  the  fatigue 
relations  to  the  nature  of  the  strength-controlling  flaw  provides  a 
strong  case  for  the  use  of  indentation  testing  in  the  evaluation  of 
crack  velocity  parameters.  Any  uncertainty  as  to  whether  the 
cracks  are  subject  to  residual  stress  effects  and  whether  the  cracks 
have  essential  point  or  line  (or  intermediate)  geometry  is  then 
eliminated.  Of  course,  such  elements  of  uncertainty  can  hardly  be 
avoided  in  ceramic  components  which  arc  to  be  placed  in  serv  ice, 
so  any  attempt  to  apply  the  above  procedure  in  reverse,  i.c.  to  use 
macroscopically  determined  crack  velocity  parameters  to  predict 
lifetime  characteristics,  needs  to  be  treated  with  extreme  caution. 
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Indeed,  the  present  resists  would  appear  to  strengthen  the  case  for 
fatigue  testing  on  the  .Mual  surface  finish  to  be  used  in  service. 

Thus  t..r  in  the  discussion  we  have  effectively  been  considering 
the  fatiguv  response  01  a  given  material/environment  system  at 
fixed  indentation  loading  conditions.  How  might  the  formalism 
developed  here  be  re.  :  ranged  to  accommodate  a  test  program 
which  ca;!‘  for  the  rem  va!  of  these  restrictions?  Ideally,  it  would 
be  convc.i'cnt  to  be  at-  to  devise  a  scheme  whereby  all  data  foi 
a  given  system  could  Is  plotted  onto  a  universal  curve,  regardless 
of  load,  such  that  comparative  evaluations  could  be  readi'y  made 
for  different  materials  Alternatively,  one  could  use  load  as  a 
control  variable  for  inscstigating  the  effect  of  flaw  size  on  the 
validity  of  crack  growth  laws.'*  Accordingly,  insetting  F.q.  (30) 
into  the  general  fatigue  relation  t,af =A’  gives’ 

t’/\\  „-(  I /v„)< (t,  i jf)"  (32) 

The  load  P.  may  now  iv  introduced  explicitly  into  (he  analysis  via 
F-q  (3)  by  writing 

(33n) 

r„=T|/f'  03b) 

where  s  t'(r,  x,,  il>..  K. )  and  tj-  rj(r,  y,,  K.  )  are  experimentally 
measurable  inert-strength  constants. 11  Then  Eq.  (32)  becomes 

i,/P.  d.Vr'i)/r„)/ltr(fll'T  (34) 

so  plots  of  t,/P?'  vs  o,Plr  in  logarithmic  coordinates  shc  .,ld  pro¬ 
duce  universal  straight  lines  for  individual  material/indenter  sys¬ 
tems  These  plots  are.  of  course,  nothing  more  than  generalized 
versions  of  those  illustrated  in  Fig.  6.  The  advantage  of  this 
scheme  is  that  it  pro- ides  the  basis  for  constructing  “master 
diagrams  in  which  the  relative  fracture  properties  of  different 
materials  ire  immediately  apparent.  Thus  the  slope  of  any  such  plot 
in  the  fatigue  region  gives  a  direct  measure  of  the  “susceptibility” 
to  delayed  failure,  in  the  manner  already  discussed  in  relation  to 
Fig.  6:  tlie  inert-strength  cutoff  likewise  gives  a  measure  of  the 
intrinsic  “toughness,"  as  reflected  by  the  parameter  ( in  Eq.  (33 a). 
Clearly,  the  materials  with  superior  strength  characteristics  will 
be  those  which,  for  a  specified  lifetime  domain,  lie  to  the  right  of 
the  diagiam.  l-urthet  details  of  this  proposed  scheme  will  be  ex- 
ploted  elsewhere  l‘* 

There  arc  in  the  analysis  several  implicit  assumptions  which 
have  not  been  given  close  attention  in  the  body  of  the  text.  These  „ 
include:  lai  that  the  solutions  of  the  fatigue  differential  equations 
are  insensitive  to  initial  conditions,  ib)  that  the  multiple-region 
effects  in  'he  crack  velocity  function  arc  negligible.  (<•)  that  the  test 
material  >>  free  of  preexisting  surface  stresses.  Reference  II  con¬ 
tains  a  detailed  discussion  of  these  points. 

Finally,  although  our  attention  has  focused  on  indentation 
ctacks  the  basic  stress  intensity  formulation  of  Eq.  (I)  might  be 
expected  to  cover  a  far  broader  range  of  flaw  types.  Similarities  in 
the  locai  residual  stress  fields  about  microstructurul  defects  in 
ceramics  (due  to  strain  incompatibilities  at  grain  or  inclusion 
boundaries)  and  indent  '.(tons  have  already  been  noted  by  Green.1' 
Extension  to  these  oiher  cases  is  accordingly  a  simple  matter  of 
reinterpreting  the  phy  deal  meaning  of  the  quantity  y,Pr  in  Eq.  ( 15) 
in  terms  of  characteristic  pressure  and  radius  parameters  which 
define  the  intensity  ami  extent  of  the  field."  Then,  provided  of 
course  that  the  reactive  chemical  species  has  access  to  the  defect 
centers  to  cause  fatigue  in  the  first  place,  the  conventional  analysis 
of  lifetime  properties  o  subject  to  precisely  the  same  modifications 
as  deser.ned  lor  contact  flaws  iri  Sections  II  and  III 


‘li  rru;  U  rurted  in  pitwnu:  th^i  it  »e  *crc  n»  t  a  dimcnMnnless  velocity 
v  *  Vr„/i,  ail ttl  Mievi  intensity  factor  K -  luj  t  »2>  uhiM  be  inverted  to  jtive 

»  * \,JCm  In  this  interpretation  the  pMttnp  scheme  to  be  [imposed  may  he  regarded  as 
rcpreventaii  .e  an  inverse  truck  velocity  function 


Failure  Strata,  o,/MPa 

Fig.  6.  Plot  comparing  static  and  dynamic  fatigue  re¬ 
sponse  for  soda-lime  glass  in  water  using  indentation  flaws. 
Data  from  Refs.  9  and  10.  standard  deviation  error  bars. 


APPENDIX 

Lifetime  Integral  and  its  Approximations 

In  this  appendix  the  assumptions  made  in  the  derivation  of  the 
lifetime  formulation  for  cracks  in  combined  residual  and  applied 
stress  fields  are  examined.  Specific  attention  is  directed  to  the 
dimensionless  quantity  in  Eq.  (26). 

AM2(r+ir,/r"]j’<V  '(I "  'rff  (A-l) 

where  the  integral  is  the  generalized  version  of  Eq.  (20)(i.e.  with 
variable  limits  of  integration).  Consideration  of  this  one  expression 
is  sufficient  to  cover  all  the  approximations  referred  to  in  Sec¬ 
tions  1H2B)  and  111.  The  treatment  here  is  taken  in  two  parts:  in  the 
first.  A,’  is  expressed  in  terms  of  the  beta  function  and  two  incom¬ 
plete  beta  functions,  and  these  functions  are  evaluated  to  determine 
their  relative  importance;  in  the  second  part,  the  beta  function, 
which  is  confirmed  to  be  the  dominant  term  in  A,’,  is  expanded  in 
an  asymptotic  series  for  large  values  of  n. 

An  important  point  to  keep  in  mind  here  is  that  accuracy  in  A’, 
although  clearly  desirable,  is  not  nearly  as  critical  as  it  is  in  the 
exponent  n'  in  Eq.  (25).  for  which  we  have  an  exact  expression 
(Eq.  (I7|).  Generally,  accuracy  to  within  a  factor  of  two  in  A,'  is 
adequate  for  most  lifetime  and  crack  velocity  evaluations. 

(!)  Beta  Function  Expression 
The  integral  in  Eq.  (A-l)  may  be  rewritten  as  follows. 

f  f  '(l-ff  "  -'#»[  f  '(l-fr  '  '<!( 

- 1  r  'd-fr  •  'd( 

-I  r  '<!-£>"  "  dt  IA-2) 

*  it 

The  first  integral  is  the  beta  function  Bin  .  n-n’)  as  previously 
defined  in  Eq.  (20).  The  second  and  third  integrals  are  the  incom¬ 
plete  beta  functions3"  Bt,  in  .  n-n  )  and  B  ,  ,,,  rit-n'. til  (us¬ 
ing  a  transformaiion  of  variables  £—*!-£’  in  (he  latter  easel. 
Equation  (A-l)  (hen  becomes 

A,'=[2(rf  I)"  7 r"  ]  [Bril'.  n-»i’t 

-B(,ln'.n-n  )-B,,  tl,(n-n‘.  n')]  <  A— 3 1 
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Table  IA,  Initial  and  Final  Values  of  Reduced 

Stress  Intensity  Factor,  |=A,/A,  for  Selected  Values 
of  Reduced  Applied  Stress,  S.=tra/<rm 

Flaw  type 

f, 

6 

Point  (r=  3)  0.9 

0.559 

0.891 

0.7 

0.390 

0.959 

0.5 

0.261 

0.986 

0.3 

0.150 

0.997 

0.1 

0.048 

1.000 

Line(r=l)  0.9 

0.342 

0.658 

0.7 

0.226 

0.774 

0  5 

0.146 

0.854 

0.3 

0.082 

0.918 

0.1 

0.026 

0.974 

Table  IIA.  Comparison  of  A,'  Evaluations 
from  Equation  (A-J)  Using  Beta  Function  Contribution 
_ Only  and  Full  Expression 


Flaw  type 

n 

ft' 

A,'  (beta) 

S„ 

a; 

Point  (r=3) 

10 

8 

1.110 

0.9 

0.780 

0.7 

1.051 

0.3 

l.l  10 

18 

14 

0.755 

0.9 

0.659 

0.7 

0.751 

0.5 

0.755 

38 

29 

0.492 

0.9 

0.482 

0.7 

0.492 

Line  (r=  1 ) 

10 

6 

2.032 

0.9 

1.181 

0.7 

1.768 

0.1 

2.032 

18 

10 

1.348 

0.9 

1.045 

0.7 

1.315 

0.3 

1.348 

38 

20 

0.864 

0.9 

0.813 

0.7 

0.864 

The  beta  function  in  Eq.  (A-3)  is  readily  computed  from  the 
more  familiar  gamma  function  to  which  it  is  related20;  generally, 

B(n',n-n')=r(n')r(n-»t’)/r(/i)  (A-4a) 

or.  more  simply,  when  both  n  and  n‘  are  integers, 

Bin'.n-n')=(n'-  1  >!(«— n'  —  I  )!/(n—  I )!  (A-46) 

Computation  of  the  incomplete  beta  functions,  however,  re¬ 
quires  prior  specification  of  the  reduced  quantities  6  and  In  the 
same  way  as  the  limiting  values  of  these  two  quantities  in  Eq.  (18) 
are  derived  from  the  defining  relation  Eq.  (15),  we  may  write, 

f=[/-/(r+l)jS0C’  (A-5) 

where  we  have  invoked  the  normalizing  scheme  of  Eq.  (23).  For 
present  purposes  we  can  adequately  illustrate  the  approximations 
involved  by  considering  the  relatively  simple  case  of  static  fatigue, 
i.e.  5„=constant;  values  of  £  and  (/  appropriate  to  initial  and  final 
crack  sizes  are  then  obtained  from  the  roots  of  Eq.(l)  at  K=Kr 
(again  expressed  in  normalized  form).  Table  IA  shows  values  for 
several  reduced  applied  stresses  for  both  point  and  line  Paws. 
These  values  can  clearly  differ  substantially  from  the  limits  of 
£=0  and  6=  I  of  Eq.  (18),  particularly  at  stress  levels  approaching 
the  inert  strength  5.  =  I . 

In  our  study  the  incomplete  beta  functions  have  been  evaluated 
using  an  adaptive  Simpson’s  rule  computer  code  for  the  integration 
limits  of  Table  IA,  although  algorithms  arc  available  for  analytical 
(but  tedious)  computation. 

Table  IIA  summarizes  the  results  of  the  calculations,  for  selected 
values  of  n  and  for  point  and  line  flaws.  The  tabulation  compares 
evaluations  of  A',  made  with  the  incomplete  beta  functions  omitted 


Table  IIIA.  Comparison  of  A,'  Evaluations  from  Equation 
(A-9)  Using  Lead  Term  and  Full  Beta  Function  Expression 


Raw  type 

n 

ft' 

[Kw/tr*  |  >*'!'  - 

V  (beta) 

Point  (r=3) 

10 

8 

0.88ri 

1.110 

18 

14 

0  670 

0.755 

38 

29 

0.46? 

0  492 

78 

1  59 

0.326 

0  335 

158 

119 

0.230 

0.233 

Line  (r=  1) 

10 

/  6 

1.447 

2.032 

18 

/  10 

1.121 

1  348 

38 

/  20 

0.793 

0.864 

78 

/  40 

0.560 

0.584 

158 

1  80 

0.396 

0.405 

and  included.  It  is  clear  that  the  terms  containing  these  incomplete 
beta  functions  make  a  significant  contribution  only  at  low  n  and 
high  S„.  Even  in  these  worst  cases  the  discrepancy  is  generally 
likely  to  be  less  than  the  factor  of  two  tolerance  level  quoted 
earlier.  Thus  in  the  region  of  larger  n  (which  appears  to  pertain  to 
most  practical  ceramics)  and  long  lifetimes.  A,'  effectively  be¬ 
comes  an  invariant  quantity,  independent  of  applied  stress. 

(2)  Asymptotic  Expansion  of  Beta  (Gamma)  Function 
In  the  spirit  of  the  approximations  above  in  which  only  the  beta 
function  is  retained  in  Eq.  (A-3)  we  have 

A,' =[2(r+ 1  )*-  7r”  ]IV)r(n-n')/r(n)  (A-6) 

The  gamma  functions  in  this  expression  may  now  be  expanded 
using  Stirling's  formula20;  e.g.  for  the  first  gamma  function 

In  r(n')~  In  (2t r/n')l,2+n'  In  n'-n'+S(n')  (A— 7) 

where  S(n')  is  the  series 

[Biil2k(2k-\j\/n'™~'  (A-8) 

with  By,  the  Bernoulli  numbers  (Bj=vi.  B4= -</*>,  etc.).20  After 
some  manipulation.  Eq.  (A-6)  reduces  to 

In  A,' —In  [8ir/(r+  l)n']1'2 

+n'  In  [«'/(«' ~2)]-(n- 16)  In  [n/(n-2)] 

+S(n')-S(n)+S(n—n')  (A-9) 

In  this  expression  the  second  and  third  terms  tend  to  cancel  and  the 
series  5  terms  become  small  as  n  increases,  leaving  the  lead  term 
as  the  dominant  quantity.  The  accuracy  with  which  this  lead  term 
can  be  used  to  represent  A,’  may  be  gauged  from  Table  IIIA. 
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